arXiv: 1509.02469vl [math.AP] 8 Sep 2015 


Michaela A. C. Vollmer 


Critical points and bifurcations of the 
three-dimensional Onsager model for 
liquid crystals 


Received: date / Accepted: date 


Abstract We study the bifurcation diagram of the Onsager free-energy func¬ 
tional for liquid crystals with orientation parameter on the sphere in three 
dimensions. In particular, we concentrate on a general class of two-body 
interaction potentials including the Onsager kernel. The problem is refor¬ 
mulated as a non-linear eigenvalue problem for the kernel operator, and a 
general method to find the corresponding eigenvalues and eigenfunctions is 
presented. Our main tools for this analysis are spherical harmonics and a 
special algorithm for computing expansions of products of spherical harmon¬ 
ics in terms of spherical harmonics. We find an explicit expression for the 
set of all bifurcation points. Using a Lyapunov-Schmidt reduction, we derive 
a bifurcation equation depending on five state variables. The dimension of 
this state space is further reduced to two dimensions by using the rotational 
symmetry of the problem and the invariant theory of groups. On the basis 
of these results, we show that the first bifurcation occurring in the case of 
the Onsager interaction potential is a transcritical bifurcation and that the 
corresponding solution is uniaxial. In addition, we prove some global proper¬ 
ties of the bifurcation diagram such as the fact that the trivial solution is the 
unique local minimiser for high temperatures, that it is not a local minimiser 
if the temperature is low, the boundedness of all equilibria of the functional 
and that the bifurcation branches are either unbounded or that they meet 
another bifurcation branch. 
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1 Introduction 

The isotropic-to-nematic phase transition of rod-like molecules is one of the 
most studied phenomena in the theory of liquid crystals. It is characterised 
by an onset of orientational order due to an excluded volume effect when 
either the concentration of molecules is increased or the temperature of the 
system is decreased. The first model describing such a phase transition was 
introduced by Onsager in 1949 m- Let p : —)> M be a probability density 
function characterising the orientation of the molecules, that is 

p{p) > 0 for all p and / p{p) dp = 1. (1) 

Using a second virial approximation, Onsager derived the tree-energy func¬ 
tional 

^{p) ■= (^kBTpip)\n{p{p)) + ^Uip){p)pip)^ dp (2) 

where U : L^(S^) —>■ L°°{S'^) denotes the interaction operator given by 

U{p){p):=[ K{p,q)p{q) dq. (3) 

The interaction kernel iL(-, •) : x —>■ M describes the excluded volume 
effect given by the interaction of one test rod with all other polymer rods in 
the system. In particular Onsager was interested in the kernel 

Ko{p,q) = \px q\ = Vi - (p-q)^ (4) 

which has been named after him. However, due to the non-analytic nature of 
this particular interaction potential, the equilibria of the energy functional 
in m cannot be computed explicitly. Subsequently, many theories and vari¬ 
ations of Onsager’s model have been formulated. One of the most popular is 
the so called Maier-Saupe potential [22] 

KMs{p,q) = ^-{p-qf- (5) 

It is based on an approximation of the mean field approach used by On¬ 
sager and assumes that the interactions between the rods are represented by 
the second moment of the probability density function for the orientation 
of the rods. A third well-known intermolecular potential is the asymmetric 
interaction potential which is also called the dipolar potential 


Kdip,q) = -p-q- (6) 

The equilibrium states of the free-energy functional above equipped with 
the Maier-Saupe kernel in have been studied extensively in the past. 
Among the first researchers who became interested in this problem was 
Freiser [12] , who extended the Maier-Saupe interaction potential to the case 
of asymmetric molecules. He proved both the existence of a first-order phase 
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transition to a nematic state and a second-order phase transition to a biaxial 
state. Constantin, Kevrekidis and Titi [3] established the equivalence of the 
equilibrium states of the Onsager functional with the Maier-Saupe potential 
in three dimensions to the solutions of a transcendental matrix equation, and 
thus derived the high concentration asymptotics in terms of the eigenvalues 
of this particular matrix. A complete classification of the equilibrium states 
in three dimensions has been provided by Fatkullin and Slastikov m and in¬ 
dependently by Liu, Zhang and Zhang [50]. Using the properties of spherical 
harmonics, both groups derived explicit formulae for all critical points and 
proved their axial symmetry. Similarly, Fatkullin and Slastikov also obtained 
results for the dipolar interaction potential in ([S|). 

One of the first approaches to the original problem involving the On¬ 
sager kernel in (|3]) has been undertaken by Kayser and Raveche m- By 
reformulating the problem as an eigenvalue problem, they derive an iterative 
scheme that allows them to compute all axially symmetric equilibria of the 
functional. Using the Fourier coefficients of the Onsager potential in two di¬ 
mensions, Wang and Zhou mm showed that there exist in fact infinitely 
many bifurcation branches with different symmetries. Revisiting the same 
question, Chen, Li and Wang |B] also proved that all solutions are axially 
symmetric in 2D. For a class of interaction potentials involving only a finite 
number of Fourier modes, Lucia and Vukadinovic HU verified the existence of 
continuous branches in two dimensions and they characterised the structure 
of the bifurcation diagram in terms of the size of the spectral gaps of the 
interaction operator. A very recent approach to the same problem has also 
been undertaken by Niksirat and Yu [23] who obtained the local bifurcation 
structure of all equilibrium states in two dimensions and the uniqueness of 
the trivial solution for high temperatures. 

However, the problem of classifying all critical points of the Onsager free- 
energy functional with the Onsager interaction kernel in three dimensions 
has not yet been addressed. In this paper, we reformulate the problem as 
an eigenvalue problem and derive a method that allows us to compute the 
eigenvalues of general interaction kernels. On the basis of this result, we 
obtain a complete set of all bifurcation points and we find the local bifur¬ 
cation structure of the first of these points corresponding to the first phase 
transition occurring at the highest temperature. In particular, we prove the 
existence of a transcritical bifurcation and we show that all critical points of 
this bifurcation branch are axially symmetric. These results about the local 
bifurcation structure of the Onsager free-energy functional are summarised 
in the following theorem. 


Theorem 1 (Local bifurcations of the Onsager free-energy func¬ 
tional) Let the Onsager free-energy functional T in m be equipped with the 
Onsager kernel Kq{ p, q) = — (p ■ q)^. Then a transcritical bifurcation 

from the trivial solution p = -^ occurs locally at the temperature 


TT 

32fcB 
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and it is uniaxial. Moreover, all other local bifurcations from the trivial so¬ 
lution occur at the temperatures 

^ _ r{s/2 + ^)r{s/2 - 1/2) 

8kBr{s/2 + l)r{s/2 + 2) 

where s G 2N and the trivial solution p = ^ is a local minimiser for all 
T > and it is not a local minimiser for r < • 

The proof of this theorem consists of two crucial steps. The derivation of 
an explicit expression of all eigenvalues of the interaction operator U corre¬ 
sponding to the Onsager kernel and the exploitation of the symmetry proper¬ 
ties inherent to the problem. In particular, the symmetry can be expressed in 
terms of a group action acting on the state space of solutions. Because we are 
only interested in solutions up to rotational symmetry, the state space of so¬ 
lutions can be reduced to the orbit space of this group action. Our approach 
involving invariant theory for groups is generally applicable to bifurcation 
problems for functions defined on the sphere and is of interest on its own (for 
details see Sections 0] and O . 

We would like to point out that even though Theorem [1] only refers to the 
Onsager kernel, our methods apply to a very general class of intermolecular 
potentials. One critical property of most physical intermolecular potentials 
is rotational symmetry. 

Definition 2. An interaction potential is rotationally symmetric, if 

K{p,q) = K{Rq,Rp) for allp,q G and for all R G SO{3). 

An important consequence of this property is the following proposition 
which is based on the basic representation theorem of simultaneous invariants 
of vectors due to Cauchy. 

Proposition 3 |291 page 29]. If the interaction potential is rotationally sym¬ 
metric, then it can be written as a function of a scalar K{p, q) = k(p-q) where 
fc : M —>■ M. 

Proof. Omitted. □ 

Remark 4. Within this paper, we use the notations K(p,q) and k{jp ■ q) 
interchangeably to denote the interaction kernel of interest. 

The following assumption characterises a large class of interaction poten¬ 
tials to which our methods can be applied to. 

Assumption 5. The interaction potential K{-, •) : x —>■ R satisfies 

(a) is continuous in both variables; 

(b) is symmetric, that is K(p,q) = K{q,p) for allp,q G S^; 

(c) is rotationally symmetric, see Definition^ 

(d) Each spherical harmonic is an eigenfunction of the interaction 

operator U in that is = XiYf^, with eigenvalue Xi so that 

I Xi I — y 0 as I —y oo . 






Critical points and bifurcations of the three-dimensional Onsager model 


5 


For a brief introduction to spherical harmonics see Appendix [A) 


In this article we restrict our attention to bifurcations of the Euler- 
Lagrange equation of the Onsager free-energy functional p = -^ equipped 
with the Onsager kernel. However, our approach is generally applicable as 
long as sufficient information of the eigenvalues A/ is available. On this basis 
our methods allow a derivation of an expansion of the bifurcation equation, 
see Remark[Tn]for more details. Applying a similar dimension reduction, one 
is faced with a similar recognition problem that needs to be solved on a case 
to case basis. All in all, a characterisation of the local bifurcation can be 
achieved using the methods presented in this research article. 

Remark 6. Similar results about the local bifurcation structure of the On¬ 
sager free-energy functional with interaction potential satisfying Assumptions 
[5](a)-(c) can also be found in the unpublished thesis of Jakob Wachsmuth EHl, 
which was only drawn to the attention of the author after completion of this 
work. The major differences between his and our work are as follows: 

— Instead of working with the set of spherical harmonics, Wachsmuth uses 
an equivalent description of this set of functions, namely homogeneous 
polynomials restricted to the sphere. This allows him to perform an elegant 
dimension reduction to a state space of two dimensions. 

— Wachsmuth does not derive an explicit expression for the set of eigen¬ 
values of the interaction operator U. He shows that the Laplace-Beltrami 
operator and U commute which allows him to conclude that the spherical 
harmonics are indeed the eigenvectors of U as well. In contrast, we derive 
a method that allows us to compute the eigenvalues of U equipped with 
any kernel admitting an appropriate Taylor expansion, see Theorem [H 

— Wachsmuth imposes the assumption that the absolute values of the eigen¬ 
values form a decreasing sequence. Assumption [5] (d) is more general and 
allows the sequence of eigenvalues to fluctuate if it converges to zero. 

— Crucially, we prove local uniaxiality of the solutions occurring at the bi¬ 
furcation point T*, while Wachsmuth only shows that the solution is in¬ 
finitesimally uniaxial. 

Last but not least, we also prove properties of the global bifurcation 
diagram which are summarised in the following theorem. 


Theorem 7 (Global bifurcations of the Onsager free-energy fnnc- 
tional) Let K{-, •) : x —>■ M 6e an interaction kernel satisfying Assump- 
tzon [S] (a)-(c). Then the uniform distribution po{p) = ^ is a critical point 
for all temperature values t. If t is such that 


SttM fl6M\ 


< T 


where M = maxp ^^§2 A'(p, g), then the uniform distribution is in fact the 
unique solution of the Onsager free-energy functional in three dimensions. 
Moreover, any non-trivial solution of the Onsager free-energy functional is 
bounded and all bifurcation branches either meet infinity or they meet another 
bifurcation branch. 
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This work is structured as follows. In Section [2] we derive a method that 
allows us to compute all eigenvalues and eigenfunctions of interaction oper¬ 
ators of the form given in ([31) which correspond to interaction kernels satis¬ 
fying Assumption [5] Carrying out a Lyapunov-Schmidt reduction, we reduce 
the infinite-dimensional problem to a five-dimensional bifurcation equation in 
Section [3] In Section |3] we reduce the dimension of the problem further by ex¬ 
ploiting the rotational symmetry of our setting which results in a simplified 
two-dimensional bifurcation equation. We solve the corresponding recogni¬ 
tion problem in Section [5] Following these results, we restrict our attention 
to the case of uniaxial solutions in Section [B] In particular, we prove that a 
transcritical bifurcation also occurs in case of the restricted problem and we 
therefore deduce that the solutions of the full problem must be uniaxial as 
well. Section [3 is devoted to a brief analysis of some global properties of the 
bifurcation diagram, such as the boundedness of all solutions, the uniqueness 
of the trivial solution as local minimiser for high temperatures and the con¬ 
tinuity of all bifurcation branches. We complete this work by summarising 
our results in Section [8] 


2 A complete description of all bifurcation points of the Onsager 
free-energy functional 

The novel contribution of this section is the full characterisation of all bi¬ 
furcation points of the Euler-Lagrange equation of the Onsager free-energy 
functional. In particular, we use the Taylor expansion of the interaction ker¬ 
nel in three dimensions. 

The Euler-Lagrange equation of the free-energy functional in ([3]) is given 


by 



where A := /c_bt. Its derivation is well known and can for example be found in 
m or [20] . The constant c is obtained through rearranging the equation and 
imposing the constraint that the orientation distribution function p integrates 
to one, see the conditions in ©• In particular, c is given by 



c = — A In Z where Z = 


denotes the partition function. Introducing the thermodynamic potential 4> : 
§2 



it follows that p{p) = Z ^ exp(—(/)(p)) and we can rewrite the Euler-Lagrange 
equation as 
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for 


Z((/)) = f exp(-(;i(q)) dq. 


The addition or subtraction of a constant to the free-energy functional in 
@ does not change its minimisers. Therefore we assume without loss of 
generality that 



k{p ■ q) dq = 0. 


Each minimiser of the Onsager free-energy functional must be a solution of 
the Euler-Lagrange equation. In particular, we prove in Proposition 1411 that 
all critical points are elements of the space The ground state of 

the Euler-Lagrange equation is given by <j)o{p) = 0 which corresponds to the 
uniform probability distribution pq{p) = ^ and solves the Euler-Lagrange 
equation for all temperatures A = fcsr. It represents the isotropic phase and 
we are interested in analysing whether other solutions exist. Mathematically 
speaking, we are interested in those values of A for which new solutions 
emerge. Locally, these are the values of A for which the implicit function 
theorem is not applicable to the operator E : (S^) —>■ (S^) given by 


:= \cl){p) - f k{p-q)exp{-(l}{q)) dq. (7) 

Z[(p) 

More details about the space (S^) and the differentiability of E and the 
regularity of critical points can be found in Appendix 1C.II In particular, E 
is infinitely many times Frechet differentiable as a mapping from (S^) —>■ 
E[^ (S^), see Lemma SH The implicit function theorem is not applicable if 
Cx : El"^ (§^) X R —>• H'^ (S^) , given by 


C-x{c^) 


dE{(j) + er], A) 
de 


= Mip) 


e=0 


47r 


dq, 


is not invertible. Hence we are interested in all non-zero solutions (</), A) of 


>^Hp) + [ Hp- Q)4‘iQ) dq = 0 

47r Jg2 

and therefore in the nullspace of the operator C\. Rearranging the equa¬ 
tion L\r]{p) = 0, it becomes apparent that this is in fact equivalent to the 
eigenvalue problem for the interaction operator in ([3]) 



k{p- q)r]{q) dq. 


In particular, we observe that for any eigenvalue /r of 17, a solution to the 
Euler-Lagrange equation is given by the corresponding eigenvector of U and 


A = - 


P_ 

Ek 


( 8 ) 
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Even though this eigenvalue problem involves only a linear integral operator, 
it is not trivial. In the following theorem, we derive an explicit expression 
for all eigenvalues of U, and thus we find the set of all possible bifurcation 
points of the Onsager free-energy functional locally around (j)o ■ In particular 
we make use of the fact that functions in (S^) can be expanded in terms of 
spherical harmonics, see Appendix E) In case of the function space (S^) 
these expansions converge uniformly. 

Theorem 8. Let x —>■ M 6e an interaction kernel satisfying 

Assumptions ^ {a)-{c) and assume that admits a Taylor expansion 

such that K(j>, q) = ar{p-qY converges for all {p, q) € {p, 9 G < 

1} and Or satisfy the condition that 


OO 


E E 

1^0 r=Z,Z+2,Z+4,... 


47r(41 + 

2(r-0/2(i(r-;))!(;+ r + 1)!! 


(9) 


Then the eigenfunctions of the corresponding interaction operator 



k{p ■ q)v{q) dq 


are given by the spherical harmonics ^ C. The corresponding eigen- 

values are given by 

E 47ras-|_2r('5 2r)! 

2’'r!(2s + 2r+l)!! 

r—0 

where s G N. 

Remark 9. Observe that we have not used Assumption [5] (d) in the state¬ 
ment of Theorem [51 Instead, we provide a sufficient condition implying As- 
sumption^{<l) by assuming that K{p,q) can be written as a convergent Tay¬ 
lor series and that its coefficients satisfy condition ^ (if the sum of the 
absolute values of all eigenvalues in is finite, the sequence of eigenvalues 
converges to zero). 


Proof. By Pi{x) we denote the associated Legendre polynomials of order I 
(for more details see Appendix IXll. Lemma 1551 in Appendix iBl states that 


= E 


{2l + l)r\ 


:Piix) 


( 11 ) 


2(^-0/2( i(r-/))!(;+ r+ 1)!! 
where the sum is taken over all r < 1 such that l—r = 0 mod 2. By assumption 

OO 

k{p-q) = ^ar{p-qY 


r—Q 


denotes the Taylor expansion of an interaction potential for all (p, q) G {p, q € 
S^IIp • 9 | < 1}. Hence an application of (ITTl) yields 


k{p-q) = Y. Y. 


{21 + l)r! 


r—0 l—r,r—2,. 


'2ir-i)n{\{r-l))\{l + r+l)\\ 


Piip-q)- 
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Using the addition theorem for Legendre polynomials [36l page 395] 


Pi{x ■ x) 


47r 

21 + 1 


Yrix)Yrix), 

m——l 


we obtain that 


k{p-q)=Y, Y. Y 


dTrOr-r! 


r=0 l—r,r—2,. 


2(-0/2(i(r-l))!(l + r + l)!! 


Yr{p)Y*Ti.<l)- 


( 12 ) 


Integrating the interaction kernel against an arbitrary spherical harmonic 
Yg{q), swapping the order of the sum and the integral and using the orthog¬ 
onality of spherical harmonics, it follows that 


n OO t n 

/ k{p ■ q)Y:{q) dq = J2 Y Y / Y^{q)Yrdq 

r=0l=r,r-2,...n.=-l 

OO I 

= Y Y E (13) 

T=0 l—r,r—2,... m— — l 

Notice that the interchange of the infinite sum and the integral needs of 
course to be verified. Viewing the infinite sum as an integral with respect to 
the counting measure, Fubini’s theorem applies if 


E 

r=0 • 


'S2 


E E c\Yr{p)Yr{q)Y:{q) 


l—r.r — 2,. 


I— — ! 


dq < OO. 


(14) 


In particular, using the Cauchy-Schwarz inequality and the fact that all 
spherical harmonics have unit mass with respect to the L^-norm, we obtain 


E 

r=0 * 


'S2 


E E c!lYr{p)Y^T{<l)Y:{q)t\p.,\<i dq 

l—r.,r—2,... m— — l 

OO I « 

^E E E \^'i\\Yr^)\ l\Yrmrm\P-<i\<A 

r=0 l=r,r-2,...m=-l 

OO I 

^Y Y Y 

r—O l—r,r — 2,... m— — l 


Moreover, it follows from [Ml Proposition 7.0.1] that any spherical harmonic 
V;™ is bounded by 


21 + 1 
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Hence, 



E E crrr(p)>^T(9)n"(9)i|p.,i<i 

l—r,r — 2,... m— — l 

E ('i' + i)="icn. 

r—0 l—r,r—2,... 


dq 


(15) 


Again we can view the two sums as integrals with respect to the counting 
measure and we may swap the order of summation according to Fubini’s 
theorem if the absolute value of the underlying function is integrable with 
respect to any order of integration. Therefore the validity of swapping the 
order of summation in m is simultaneously proved when we establish our 
original claim in (O. Therefore a basic condition that suffices to be proved 
in order to guarantee the validity of the interchange of the integrals in m 
is given by 


E E (4Z + l)^/^|cri<oo. (16) 

1=0 r=l,l+2,l+4,... 


In particular, one needs to show on a case to case basis that the coefficients 
|c[| decay faster than In Lemma [551 we prove this condition in case of 
the Onsager kernel. For all other cases this condition is assumed to hold (see 
the statement of the theorem). 

Having established m, we deduce that 


kip-q)Y:{q)dq = Y, 


r—0 


47ra,+2r(5 + 2r)! 
2’-r!(2s + 2r + l)!! " 


(17) 


From the last equation we may conclude that the eigenfunctions of the inter¬ 
action operator Uri{p) = Jga k{p-q)r]{q) dq associated to the class of two-body 
interaction potentials k(p ■ q) described in the statement of this theorem are 
given by the spherical harmonics and that their corresponding eigenvalues 
are given by 


E 47ras+2r(s + 2r)! 
2^r\{2s + 2r+l)\\ 


with s G N. 


□ 


Corollary 10. The eigenfunctions and eigenvalues of the interaction oper¬ 
ator U associated to the Onsager kernel ko{p ■ q) = -^/l — {p ■ qY are given 
by the spherical harmonics 

{y” G L^(S^) : s G N even ,—s<n<s} 


and 


do{s) 


7i-r(s/2-i-i)r(s/2-i/2) 

2 r(s/ 2 -|-l)r(s/ 2 -|- 2 ) 


if s is even 
if s is odd 


0 
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Proof. Based on the Taylor expansion of the square root function for 

all X G (—1,1), it is easy to see that the Taylor expansion of the Onsager 
kernel is given by 


ko{p-q) = TT-^ 


(2r)! 


r=0 


(1 - 2r)(r!)2(4'-) 


ip-q? 


for all p, g G such that \p ■ q\ < 1. We deduce that 


ko{p-q) = Y^ ar{p ■ qY with Or = i (i-0(5!)^(20 

r=0 I® 


if r is even 
if r is odd 


and thus using m 

ko{p-q) 


= E E E TT 


47r(2r)!2r-(p)r*r(9) 


Eli (1 - 2r)(r!)2(4’-)2-V2(^ _ ;/2)!(Z + 2r + 1)!! 
47r(2r)!2y-(p)r*^(g) 


r=0 I—'. 
oo r 21 

= EE E 

r=0 Z=0 m— — 2l 


(1 - 2r)(r!)2(4’-)2'-'(r - l)\{2l + 2r+ 1)!!' 


(18) 


Having expressed the Onsager kernel in terms of spherical harmonics, we 
may apply m in Theorem [5] in order to get 


ko{p ■ q)YY{q)dq 


/S2 


EE E 

' r—0 /—0 m— — 2l 
oo oo 


4^(2r)!2r-(p)E*^(g)n"(g)l|p.,l 


<1 


(1 - 2r)(r!)2(4'-)2’'-'(r - l)\{2l + 2r + 1)!! 


dq 


E E (1 _ 2r)(r!)2(4’')2'-^(r - l)l{2l + 2r + 1)!! E ^2T(p)<5mn<5(2i) 


l—O r—l 
■^oo 


= -2Z 


= J ^r=s/2 (l-2C(r!)^(4C2-l°V(r-^/2)!(s+2r+l)!! ^"(P) even 

1 0 if s is odd 


where we used the orthogonality of the spherical harmonics and Lemma [35] 
in which we verify the condition in (1161) that allows an interchange of the 
sum and the integral. 

In fact, using the theory of generalised hypergeometric functions, it can 
be shown that this expression for the eigenvalues can be written in a closed 
form (for details see Lemma 1371 in Appendix iBl) . It reads 


Po{s) 


7vr{s/2+^)r{s/2-l/2) 

2 r(s/ 2 +l)r{s/ 2 + 2 ) 

0 


if s is even 
if s is odd 


□ 
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Based on Corollary [10] and ([8|) , we therefore conclude that a set of all 
possible bifurcation points is given by 



Remark 11. We will show in Theorem IMI in Section^that in fact any point 
Ao(s) is a bifurcation point. 


3 The bifurcation equation of the Onsager free-energy functional 

Having established a general method for deriving the bifurcation points of the 
interaction operator, we direct our attention to the corresponding bifurcation 
equations. We would like to mention that this procedure is quite general 
and may be applied to a very large class of interaction kernels that satisfy 
Assumptions|S](a)-(c). However, in the following analysis we will restrict our 
attention again to the Onsager potential. 

The presentation of our results is divided into two parts. We begin by 
giving a brief introduction to the Lyapunov-Schmidt decomposition, which 
is our main tool for computing the bifurcation equation, see Section [01 In 
Section 13.21 we focus on the practicalities that are involved when we apply 
the Lyapunov-Schmidt decomposition to the Euler-Lagrange operator given 
in ©. In particular, Sectionis divided into four parts: a reformulation of 
the problem in the language of Section IXTl see Section 13. 2. II an algorithmic 
procedure that allows us to fulfil the decomposition, see Sections 13.2.21 and 
13.2.41 and the presentation of an algorithm that allows the fast computa¬ 
tion of products of spherical harmonics in terms of spherical harmonics, see 
Section 13.2.31 


3.1 The theory: The Lyapunov-Schmidt decomposition 

This brief introduction to the Lyapunov-Schmidt decomposition is based on 
[5]. The main idea of the Lyapunov-Schmidt decomposition is the reduction 
of a possibly infinite-dimensional bifurcation problem to a finite-dimensional 
one. Let 


F{w,X) = 0 


(19) 


be the equation of interest with u; G AT, A G M and F : W x R —>■ R where X 
denotes a Banach space. Without loss of generality we assume that 


F(0,0) = 0. 


The linear and non-linear parts of m decompose as 
F{w, A) = A) -t- TZ{w, A), 
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where 


,, , ar(o + e,,o) 
-ai- 


€=0 


( 20 ) 


and where TZ denotes its remainder. We assume that £ is a Fredholm operator 
of index zero, that means that the kernel N(L) has finite dimension d, the 
range R{L) has finite codimension r and that both dimensions d and r agree. 

The idea of the Lyapunov-Schmidt decomposition is to reduce the di¬ 
mension of the equation and its solution by projecting it onto the kernel of C 
which is, due to the assumptions on C, finite-dimensional. Let P : X ^ 
denote the projection onto the nullspace of £ and let (1 — P) be the projec¬ 
tion onto its complement. Then m is equivalent to the system of equations 


PF{u + v,X)=0 (21a) 

(1 - P)£(w-hu,A) = 0 (21b) 

where u := Pw denotes the projection of the solution onto N{C) and v := 
(1 — P)w its projection onto the complement. The second of these equations 
can uniquely be solved for v using the implicit function theorem and the 
solution v{u, X) can then be plugged back into the first equation, which yields 
the bifurcation equation 

fiu,X) := PTZ{u + v{u,X),X). (22) 

Note that we dropped the first term appearing in (I21al) because £ is linear, 
and thus it follows that PC{u -I- f) = 0. Similarly, we can drop the linear 
part of (Eg. The solutions of (1221) are equivalent to the solutions of our 
original problem in m , see [5], because it only depends on u, which is an 
element of the finite-dimensional kernel of £, and thus reduces the possibly 
infinite-dimensional problem to a finite-dimensional one. 


3.2 Practicalities: An algorithmic procedure to derive the bifurcation 
equation 

According to the previous section, the computations for the derivation of the 
bifurcation equation in (1^^ are divided into two parts: 

Step A. An application of the implicit function theorem to (I21bl) that 
gives us an explicit expression of v in terms of u and A. 

Step B. Plugging v back into (EUl) which results in the bifurcation 
equation. 

By using the implicit function theorem in the first of these two steps, 
it is often not possible to derive an expression in closed form. Instead, we 
obtain a Taylor expansion of v{u, A) (which is justified by the implicit function 
theorem [3l Theorem 2.3]) and its coefficients are obtained by matching those 
of the same order that occur in (l21aL The resulting expression for v{u, A) can 












14 


Michaela A. C. Vollmer 


then be plugged into the expansion of TZ. An application of the projection V 
finally yields an expansion of the bifurcation equation in (1221) . 

In order to characterise the bifurcation occurring at A 2 , it is often suffi¬ 
cient to obtain an expansion of (EH) up to a particular order. We know that a 
certain order is sufficient by looking at the so-called recognition problem that 
corresponds to the problem at hand. More details on the particular recogni¬ 
tion problem that we need to apply in order to solve the bifurcation equation 
can be found in Section [SJ where we will show that it suffices to compute the 
bifurcation equation up to third order. 

In the following, we will first restate the problem in the language of Sec¬ 
tion [23] before directing our attention to the first of the two main steps men¬ 
tioned above in Section 13.2.21 In order to make all computations tractable, 
we write all expressions in terms of spherical harmonics. An essential tool for 
the comparison of coefficients with matching order is an algorithmic proce¬ 
dure that allows us the fast computation of products of spherical harmonics 
with a computer. This will be the subject of Section r3. 2. 31 We conclude this 
section by performing Step B above which yields the bifurcation equation in 
five dimensions. 


3.2.1 Reformulation of the problem 

In order to find the minimisers of the free-energy functional, we consider its 
Euler-Lagrange equation and thus the corresponding Euler-Lagrange opera¬ 
tor 


A) = {\s + X)(t){p) - [ k{p ■ q) exp(-(/)(g)) dq. (23) 

Z[9) 

Remark 12. Notice that (j23p is an update of the Euler-Lagrange operator 
in 0 which has been translated by a factor Xs in order to ensure that the 
assumption £’(0,0) = 0 holds. 

The constant A^ denotes the bifurcation point of interest. According to 
dS]), it is given by 


where fj-s denotes the eigenvalue of order s of the interaction operator as¬ 
sociated to the Onsager kernel, see Corollary [TUI for details. The variable A 
denotes the deviation from this bifurcation point and will act as bifurcation 
parameter. Having this translation of the problem in mind, we are interested 
in bifurcations around the point {(f, A) = (0,0). In particular, we are looking 
for the first phase transformation that occurs when the temperature is de¬ 
scending. Hence we consider the largest bifurcation point which is, in case of 
the Onsager kernel and according to the result of Corollary [121 given by 

M2 _ tt 
^ dTT 32 ■ 


In order to derive the decomposition of the above functional into its linear 
and non-linear part, we consider its Taylor expansion E up to fourth order 
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in both variables (j) and A around the point (0,0), namely 

:={X 2 + X)(j){p) + ^ [ k{p ■ q)(j){q) dq + [ k{p ■ q)(j){q) dq 

An Jg2 7g2 

- ^ f k(P- q)4>'^{d) dq + [ k{p ■ q)(j){q) dq 

Stt J§2 7g2 

- k{p ■ q)(i){q) dq - Hp ' <l)4>^{q) dq 

+i i/*” ■ i '''' 

OAtT 7g2 Jg2 

+L ‘''' “ i '''' 

+ I «!■ ■ *! + ^"‘seS X ''I'’ ■ 

(24) 

According to (I2()L its linear and non-linear parts therefore decompose as 

£(</)) = A 2 <(i(p) + ^ [ k{p- q)(t){q) dq 
47r 7g2 

and 

n{4>.x)-.= H<t>A)-c{4>), (25) 

respectively. The fact that £ is a Fredholm operator of index zero is a con¬ 
sequence of the following corollary in [T] . 

Corollary 13 [3 Corollary 4.47]. Let T, K : X ^ y be linear operators and 
assume in addition that T is a Fredholm operator and that K is compact. 
Then T + K is also a Fredholm operator with index(r -|- K) = index(T). 

In Lemma [40] in Appendix 1C.II we prove that the interaction operator 


U{(l))ip)= [ k{p-q)e^p{-<p{q)) dq 
JS2 

is a compact operator for all interaction kernels satisfying Assumption |S] In 
order to prove compactness, it is crucial that the operator satishes Assump¬ 
tion [5] (d). Schematically, C can be written as 

C = C Id+ K. 

Hence in order to apply Corollary [131 we only have to show that the identity 
is a Fredholm operator of index zero. Let us make the following observations 
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(a) Id is bounded and linear 

(b) A/’(Id) = {0} and coker(Id) = {0}. 

Hence the identity is a Fredholm operator of index zero and thus so is C as 
claimed. 

In view of the results of Section [51 it is easy to conclude that the kernel of 
the operator £, denoted by N (£), is the eigenspace of the spherical harmonics 
corresponding to the degree of the bifurcation point of interest, which is in 
our case A 2 . Thus, 

N{C) = {^ 2 ™ : -2 < m < 2} 

which is a five-dimensional space. It follows that the projection P onto N{C) 
is given by 


2 

P : L^{S^) where Pw := ^ (w, ^ 2 ’”)^=^^' for all w G 

m—— 2 

Let u{p) := P4){p) and v(p) := (I — P)4){p) such that u(p) -l-u(p) = w(p) and 
let us further define Um and vi^m such that 

2 00 I 

u{p) ■■= X] '^rnX^ip) and v{p) ■='^ vi^„,Y/^{p), (26) 

m— — 2 Z=0 m= — l 

respectively. We will see that it will be an essential step to write all ex¬ 
pressions in terms of spherical harmonics as this simplifies the actions of all 
operators. In particular, we may deduce from Theorem [5] that 


Uivn = [ Hp- q)Yriq) dq = piYrip). 

JS2 


Accordingly, we conclude that 

CiYn = ( a . + £) Yi 


and similarly. 


c-\Yr) = 


Att 


47rA2 -I- Pi 


-Yr 


(27) 


Since the set of all spherical harmonics is an orthonormal basis for the space 
L^(S^) and orthogonal for PP' (S^), these operators can in fact be viewed as 
multiplication operators on the space of functions L^(S^), such that Cw{p) = 
(•^2 + w{p) for all w G L^(S^). This leads to an enormous simplification 
of our subsequent computations. 
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3.2.2 Step A: An application of the Implicit Function Theorem 

It is straightforward to show that the implicit function theorem is applicable 
to 


Cv + (l-P)n(u + v,\) = Q. (ITThll 

The operator C is invertible in R{C) and hence its inverse is bounded. Since 
Zl„72.(0,0) = 0 by definition, the implicit function theorem is applicable. 
However, writing down the exact solution v{u, A) is not possible in practice. 
Instead, an algorithmic procedure is required which has been presented in 
[5] . It is mainly based on writing as a Taylor expansion in terms of both 
variables u and A which is justified by the differentiability of the implicit 
function, see implicit function theorem [3] Theorem 2.3]. In particular, we 
expand f = (I — P)(/) up to fourth order by 

v{u,X){p):= ^Vi^j{u,X){p), 

where the terms Vij denote terms of order in u and of order in A. 
Furthermore, we identify u as before with 

2 

u(.p) = Y UmYf^ip). 

m—— 2 

By definition v G N{£)^ and thus, without loss of generality, we assume that 
Vij € (iV(£))-*- for all i,j. Plugging this expression into (I21bll 

C{v) = -(I - P)TZ{u + v,X) 

and matching the terms of the right order on both sides, we obtain expressions 
for each of the Vij. We will execute this calculation for the first couple of 
terms. 

Remark 14. Because the spherical harmonics form an orthonormal basis 
for the Vij can be written as an infinite series expansion. In fact, 

when considering the first few steps of the procedure, we will see that this 
expansion is finite. 

Since there does not exist any term that is constant in both variables, we 
conclude that 

^ 0.0 = 0 . 

As a second step, we consider Di^o which only depends on u. Hence we only 
consider terms of order one in u on both sides of (I21b|) taking the expansion 
of E in (IM)) and the definition of TZ in (1^51) into account. All terms on the 
right hand side depend on A or on higher order terms of u, so they can be 
neglected when we match the coefficients. Thus, we obtain 


k{p ■ q)vi^o{q) dq = 0. 
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This is equivalent to solving /h)i,o = 0. Because we are assuming that Vij € 
iV(£)-*-, we conclude that vi^o = 0. Similarly, one can deduce that 


Do,i = 0. 

As the first non-zero term, we consider ^V 2 ,o- Matching the order of terms 
on both sides of (I21b|) . we obtain 

A If 

7:^2,o{p) + i- / k{p- q)v 2 ,oiq) dq 

= -(1 - P) ■ d)u{q) dq - j^ Hp ■ qWiq) dq^ . 

This yields an explicit expression for V 2 fi if we bear in mind that the operator 
C and its inverse can actually be interpreted as simple multiplications, 
see (Ell). However, in order to be able to apply this multiplication to the right 
hand side, we need to write it in terms of spherical harmonics. Recall that 
u itself is already given as an expansion in terms of spherical harmonics, see 

(uni), 


2 

u{p) = X] u„,Y^{p). 

m— — 2 

Hence we only have to expand the terms 


/ kip ■ q)uiq) dq 

J§2 


uiq) dq and 


Hp-qWiq) dq. 


In fact, 

/ kip ■ q)u{q) dq / uiq) dq = 0 

because the second factor, being the sum of integrals of spherical harmonics 
of degree I = 2, is zero. In case of the second term, we observe that we have 
to compute which is based on computing products of spherical harmonics. 
In order to fulfil this task, we have developed an algorithm that is presented 
in Section 13.2.31 Applying it to the case above yields an expansion of 
in terms of spherical harmonics. We are now in the position to apply the 
interaction operator U, the projection I — P and the operator C~^ which, 
according to (EZl), reduces to a simple multiplication in this case. We obtain 


V 2 ,Qip) = -C M (I - P) / kip- q)u^iq) dq) 
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-h 


-h 


64 


Uy,-Hp) + 


64 


l-lU-2Y^ ■^{p) + 


•\/l57r" 

448 


U0M-2X4 ^{p) 


r3/2 


'(p) + ^U2U-2Y^ip) + ‘^(p) 


15^3 


+ 


+ 


224 


iSttS 

2 '' 


224 

■v/TStt^ 

448 


-u^iUqY^ ^{p) 


o 3/2 3/2 

+ —u.^u^Y^ip) 


448 


112 


uoUiYlip) 


448 


W_1M2F4 (p) 


224 

— TT^ 




M0M244^(p) 


64 


-MlU2y/(p) 


14 


■^yiip)- 


64 


This concludes the computation of 1 ) 2 , 0 ■ The computation of the other terms 
follows the same pattern. However, because they are more complex and te¬ 
dious, we will omit them. In order to give a good overview of the individual 
steps for each case, we present a summary of the procedure in form of a 
general algorithm. We would like to conclude this section by remarking that 


Algorithm 1: Matching the terms in Equation (I21bll 
1: Input: Fix Vij such that 0 < i -I- J < 4. 

2: Write down ']Z{u + v) keeping only terms in u and all terms in Vij up to i’th 
order and in A up to j’th order. 

3: Expand TZ(u + v) in terms of spherical harmonics using the results of Section 

[ 3 :^ 

4: Apply the projection (1 — P) by dropping all spherical harmonics of degree 
I = 2 . 

5: Apply -jC-^ by multiplying 7^(a -b v) by 
6: Output: Vij = + v). 


the expansion of v resulting from this procedure is in fact finite and that this 
is always the case when considering a finite expansions of an equation. The 
reason for this is that rt is a linear combination of spherical harmonics that 
span the kernel of C which is itself finite. Considering products of up to four 
terms in u and Vij, we can therefore only reach spherical harmonics up to 
degree I = 8^ since the highest degree for Vi j with i + j = 4 can only be 
1 = 8 . 

Remark 15. The operator C depends on the bifurcation parameter and the 
eigenvalues of the interaction operator U. For any kernel satisfying Assump¬ 
tion 0 the expansion of v would only depend on finitely many eigenvalues 
corresponding to spherical harmonics with low order. That means the result¬ 
ing bifurcation can then be studied with the methods presented here. 
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3.2.3 Products of spherical harmonics 


The crucial step in matching the coefficients in (I21bl) is to expand products 
of spherical harmonics again in terms of spherical harmonics. It is known 
that this can be achieved using Clebsch-Gordan coefficients which arise in 
angular momentum coupling m However, it is very hard to compute them 
explicitly because it takes too much computing time. Therefore we devised 
an algorithm for the computation of products of spherical harmonics with 
short run time on the computer. The exact algorithm that we implemented 
on the computer is given in Appendix 1C.21 while its crucial ideas are subject 
of this section. 

We first recall the definition of spherical harmonics (see Appendix |2| 
Yri<p, 9) = A^i™e™®Pr(cos ^p), -I < m < I 


where p and 9 denote the polar and the azimuth angle corresponding to 
the unit vector p S S^. The functions P™ are called associated Legendre 
polynomials and their exact definition is also given in Appendix Thus, a 
product of two spherical harmonics is given by 

Yrip, 9)Y^{p, 9) = Vr(cos(p))Pp«(cos(p)). (28) 

Notice that we will abuse the notation and drop the (^, 9) or p dependence 
in due course. Due to the factor we observe that the order of each 

spherical harmonic occurring in an expansion of the product in (1281) has to 
be of order m + g while its degree may be arbitrary. On the basis of this 
observation, we would like to investigate products of associated Legendre 
polynomials. In fact, we use the following recurrence formula from |33j 

{l-m+ = {21 + l)xPr - {I + m)Pfli. (29) 


Rearranging it gives 


pr = 


21-1 


-xPi-i - 


I — 1 + m 
I — m 


P, 


1-2 


and thus. 


pmpq_‘^ _ItP"* P'3 — ^ ^ 

-O - I _ 




pm pq 
'^l-2^p ■ 


Applying the original recurrence rule in (12911 again, but this time to the 
product xPp^q, yields 


T-pg — 

P l-m 


21 — ifp — q+1 
2p+l 


pQ pm 
^p+lM-1 


2p+l 




l-rn 


Using the definition of spherical harmonics, we deduce that 

ymyq ^.l-lp-q+l Ni^rn P^P,q ym yq 

I p -l-m 2p+\ ^>+1 

I 2^—1 p + g Npm Np^q „ 

l-m2p+l ^{p-i),q 

/ - I + TO W,m 
I - TO 


^^Tn ^rq 
^l-2^p 
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and thus, we observe that we reduced the degree of the spherical harmonics 
yj™ by one or two. We repeat this procedure until we reach either a term of 
the form y7-fi (which has to happen eventually). Hitting first, we 

use another recurrence rule from |33j . namely 



so that we obtain in any case YfY^. Again we apply a recurrence rule 


P® = (-l)"(2s-l)!!(l-a:")*/2 


which yields 


r/y/ = (-i)*(2s -1)!!(1 - x^y/^Ns,sYp^. 


This last product can then be resolved by using 



repeatedly if necessary. Again this recurrence formula arises from using re¬ 
currence formulae for associated Legendre polynomials. 

The computation of the first couple of steps of the algorithmic procedure 
above illustrates how an expansion of a product of two spherical harmonics in 
terms of spherical harmonics can be obtained. The general method involves 
a couple of other cases to start with, such as spherical harmonics of negative 
order for example, which also have to be taken into account and for which 
special rules need to be defined even though the general method works exactly 
like this. A detailed list of all rules that are necessary to cover all cases that 
might occur is given in Appendix 1C.21 In order to implement this algorithm, 
it is important to use a conditioned replacement rule that makes it faster. 
That means that the algorithm itself remembers cases that it has computed 
already so that it will stop once it hits a known target. 

3.2.4 Step B: Derivation of the bifurcation equation 

Having found an explicit expression for v{u, A) in Section 13.2.21 we are now 
in the position to plug this expression into (1^ which finally gives us an 
approximation of the bifurcation equation up to fourth order 


/(m. A) := P'R(u + v{u, A), A). 


Remark 16. One can show that the bifurcation equation corresponding to 
the Onsager free-energy functional is smooth, see Lemma 1431 

Remark 17. In fact, we used a fourth order Taylor approximation of the 
Euler-Lagrange operator and a Taylor approximation ofv(u,X) up to a joint 
fourth order in both variables. Thus, we obtained an approximation of the 
bifurcation equation of at least fourth order. More precisely, the bifurcation 
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equation f = TZ{u + v{u,\),\) = 0 is approximated based on the following 
expansions of v and TZ 


v(u,A) — v(u,X) + O I ll'*^llff 2 (g 2 )|A| 

Vi=l 

A) = A) + O (Eihi'h. (S=)l^l 


5—i 


5-i 


The following calculation shows that f agrees with f up to fourth order 
f =TZ{u + v{u, A), A) 

=R{u + v{u, A), A) + O 11'“ + '^)llff 2 (s 2 )|A|® *1 


vi=l 
/' 6 


—']Z{u + v{u, A), A) + O I ^ I |u| |ff 2 (s 2 ) IA|® 


+o El 


I//2(§2)| 


\2=1 
|5 —2 \ 


^2=1 


f{u, A) =f{u, A) + O ^ 11(/.|r^ 2 (g 2 ) IA|^-* . 


^i=l 


The bifurcation equation can also be written as an expansion in spherical 
harmonics of degree 1 = 2, and therefore admits the form 


f{u,X)ip)= fmiu,X)Yf^{p). 


(30) 


m= —2 


The solutions of this equation are equivalent to the solutions of the Euler- 
Lagrange equation of the Onsager free-energy functional in ©. In order to 
find the zeroes of (|30l) . each of the coefficients for — 2 < j < 2 needs to be 
zero. Therefore we seek the solutions to a system of five polynomials, each 
depending on the six variables m_ 2 ,u_i, uq, ui,U 2 and A which reduces the 
infinite-dimensional state space of our problem to the ten-dimensional 

one i.e. C®. 

However, we are only interested in real solutions to this equation, so we 
can reduce the dimension further by restricting the equation onto the space of 
real spherical harmonics with real coefficients. In particular, the real spherical 
harmonics are defined as 

'-i-irr-n ifm<o 

Yi^ if m = 0 . 

+ if^>o 


Yl^rn •— ^ 
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Thus, if the function u{p) in (IMl) is real-valued, its coefficients 


can be written as 


u := (u_ 2 ,n-i,Uo,Ui,U 2 ) 


/ 7f 0 0 0 


u = T(a) := 


V2 


0 ^ 0 
0 0 10 0 
0 

0 0 0 


— 0 - — 
V2^ V2 


0 


\ V2 



/ 0-2 \ 


a_i 


Oo 


Oi 


\ 02 / 


(31) 


where T : —>■ C® for some a := (a_ 2 ,..., 02 ) S K®. Using this transforma¬ 
tion, we define the following real version of the bifurcation equation which 
corresponds to the set of real spherical harmonics {d 2 _m : —2 < to < 2} 


Aeai(a,A) :=r-V(Ta,A). 


Equivalently, we will use the expression 


/,eai(a,A) :=r-i/(Ta,A). (32) 

in order to denote its approximation up to fourth order. Due to its compli¬ 
cated form, we will not state the bifurcation equation in terms of real spher¬ 
ical harmonics explicitly in this section. Instead it can be found in Appendix 

ICAl 

Remark 18. From now on, we will refrain from referring to (1321) as bifur¬ 
cation equation in terms of real spherical harmonics but we will also simply 
refer to it as bifurcation equation. 


4 Reducing the dimension of the state space of the bifurcation 
equation 

The most common approach to investigate the local structure of a system 
of polynomial equations at a given solution point is to use its Groebner ba¬ 
sis. A Groebner basis is the generating set of an ideal in a polynomial ring 
over a field K[xi,X 2 , ■ ■ ■ ,Xn]. The advantage of this method is that it can 
be used in order to find the dimensionality of the solution set at a given 
point. However, due to the lengthy form of the bifurcation equation in (I32L 
this method is not directly applicable. Instead, we make use of the symme¬ 
tries of our problem in order to reduce its dimensionality. In particular we 
consider the Onsager free-energy functional with a rotationally symmetric in- 
termolecular two-body potential. That means, that if we can find a rotation 
that translates one solution of the bifurcation equation (1321) into another, we 
can consider both solutions to be equivalent. 
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In order to represent an arbitrary rotation matrix with 

respect to the Euler angles ipn, '4’R ^-nd d/j, we use the so-called y-convention. 
In this convention, a rotation R{(pR, ipR, Or) is given by the product 


R{^R^'4’R, Sr) = Rz ■ Rx ■ Ry 

where the first rotation is by an angle ipR around the z-axis, the second 
rotation is by an angle Or around the y-axis and the third rotation is by an 
angle ipR around the z-axis. In particular, the three matrices are given by 


/l 0 0 \ 

i? 3 ; = 0 cos{0r) sm{0R) ,Ry 

yO -sm{0R) cos{9r) J 

( sm{ipR) -cosipjR) 0\ 
cos(-!/'ii) sin(V’fl) 0 . 

0 0 l) 


( -sm{ipR) cos{ipr) 0 \ 

- cosi(pR) - sm{(pR) 0 , 

0 0 l) 


In contrast to the rotation of points in the state space, the rotation of a 
function u{p) written in terms of complex spherical harmonics of order I = 2 
corresponds to a linear transformation of its coefficients u = (u_ 2 , ■ • ■ ,U 2 )- 
This transformation is given by the 5x5 Wigner matrix D((Pr^'iPr,0r). In 
particular, this relationship is formulated as 



{R{ipR,tpR,9R) -p) 


V 


{D{pr, IpR, 0R)u)my2 


{P) 

(33) 


for all p G [53] . If we are given two vectors a, b G each representing 
the five coefficients of a real-valued solution written in terms of real spherical 
harmonics, then we say that a ^ b if and only if there exists a tuple of 
rotation angles {ipR,ipR,0R) such that 


M{^r,P^r, 9R)a := T-^D{PR, p,R, Or) • Ta = b (34) 

where T denotes the isomorphism given in (I3T]) . In other words, M is the real 
version corresponding to the complex Wigner matrix D{(pR, ipR, Or) given in 

(|33|. 

Having introduced these concepts, we are now in the position to establish 
symmetry properties of our bifurcation equation. 

Proposition 19. The bifurcation equation /real is equivariant with respect to 
the Wigner matrices M{ipr,'iPr,Or), that means that for all M(ipr,iPr,Or) 
and every a G M®, 

freal{M{(fR,1pR,0R)a,X) = M {(fiR, fjR, 0 r) freal{a, X). 

Proof. Let R{(Pr,'iPr,0r) be the rotation matrix that corresponds to the 
Wigner matrix M(pr, ipR, Or). For ease of notation, we drop the dependence 
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of R on the Euler angles and 9r. Defining R(j){p) := (t>{Rp), we see 

that 

E{R4>, A) =\(t){Rp) - ' [ k{p ■ q) exp(-(/)(i?g)) dq 

Z{R(p) Js2 

=X(j){Rp) - j - \ / k{p ■ q) eyip{-(j){Rq)) dq 

Jg 2 exp(-0(i?s)) ds Js 2 

=\(l){Rp) - -j - / . '' [ k{p ■ R*q) exp(-0(9)) dq 

Js 2 exp(-</>(s)) ds Jg 2 

=\(j)[Rp) - -J - / '' , / k(Rp ■ q) exp{-(l){q)) dq 

Jg 2 exp(-0(s)) ds Jg 2 

=RE{(j), A). 

Hence the operator E is equivariant with respect to the Wigner matrices. By 
definition 

/ := P7^(u^-v(u,A), A) 

where 72. denotes the non-linear part of E and is therefore also equivariant. 
The projection P maps all terms onto the spherical harmonics of degree I = 2 
and therefore also preserves the equivariance as well as the approximation 
of the equation up to fourth order. The restriction to the subspace of real 
solutions is also a linear transformation and therefore we can deduce that 
/real Is equivariant with respect to the Wigner matrices M[ipr,iPr,9r). □ 

Knowing that / is in fact equivariant with respect to rotations, we can 
restrict our attention to one representative of each class of solutions that can 
be translated into each other by a rotation. In other words we are interested 
in the orbit space that corresponds to the action of SO{3) on the space M®. 
To illustrate this idea, consider the example of SO{2) acting on The orbit 
space in this case can then be taken to be the non-negative part of the x-axis 
and is therefore one-dimensional. 

The aim of this section is to show that the orbit space of SO{3) acting on 
K® can in fact be reduced to two dimensions. In particular, using invariant 
theory for groups, we can prove that the space 

5 := {(0,0,x,0,y)|a;,j/G M} (35) 

contains at least one representative of every orbit. 

In particular, one can show that the generators of the ring of G-invariant 
polynomials for any compact Lie-group G separate the orbits; thus for any 
two distinct orbits E and E', there exists at least one of the generators taking 
different values on E and E' [H Appendix C]. In the case of the group action 
given in (13411 there exist two invariant polynomials Ii and I 2 generating the 
ring. The fact that these generators separate the orbits can be reformulated 
as follows. If 


(/i(a),/ 2 (a)) = (/i(b),/ 2 (b)). 
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then there exist rotation angles ^pn, ipu and Or such that M{ipR, ipR, = b 
and hence a ^ b. This allows us to check that S contains at least one 
representative for every orbit by verifying that 

{/i(a),J 2 (a) |aeM5}={/i(a),/2(a) |ae5}. (36) 

In other words, we show that the state space can in fact be reduced to two 
dimensions. 

In [ 30 ] Wachsmuth achieves a somehow similar two-dimensional reduction 
by the interesting device of considering the set of homogeneous polynomials 
restricted to the sphere, which is an equivalent description of the set of spher¬ 
ical harmonics. 

This section is structured as follows. In Section [4. 1 1 we find an isomorphic 
representation of the group action, called the Cartan representation, and we 
derive an explicit way to separate its orbits based on invariant polynomials. 
Using the isomorphism between the Cartan representation and the group 
action given in (IMll which is based on the Wigner matrix, see (IHTll . we obtain 
similar expressions for the invariant polynomials with respect to the group 
action in (IMl) . We conclude this section by proving that the reduced state 
space S contains a representative of every orbit, see Section IT^ 


4.1 The separation of orbits of the group action SO{3) acting on R® 

There exists a unique irreducible representation of the group action of SO{3) 
acting on [13]. Since the complex and real representations based on the 
Wigner matrix given in Appendix ID. II and ID. 21 are irreducible, it is isomor¬ 
phic to any other irreducible representation of SO{5) acting on the five¬ 
dimensional real space. One of these isomorphic representations is called the 
Cartan representation, see Appendix ID. 31 The explicit isomorphism between 
the two group actions is given in Lemma |33l in Appendix ID. 41 The Cartan 
representation is given by 

p : 30(3) X su(3)/so(3), (S', X -H so(3)) ^ SXS-^+so{3) 

where su(3) and so(3) denote the Lie-algebras associated with the Lie groups 
SU{3) and SO(3), respectively. The representation above then arises as the 
adjoint representation of so (3) [5]. In this setting the matrix X G su(3)/so(3) 
is a traceless skew-Hermitian matrix, that is X* = —X. In particular, if we 
write X = U + iV, with U and V both being real 3x3 matrices, then we 
can make the following observation 

X* = U* + {iVy = U* - iVK 

Thus, in order for X to be skew-Hermitian, U needs to be a traceless skew- 
symmetric matrix, thus U € so (3), and V needs to be traceless and symmet¬ 
ric. It follows that 

su(3)/so(3) = {X-hso(3)|X gsu(3)} 

= {iV -f so(3)|tr(U) = 0 and U‘ = V}. 
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Let us choose the following basis for this space 


E = 



1 0 0 \ 

0 

0 

0 

00 0 

, 62 = 0 1 0 

0 0-1/ 

yo 0 -1 


/o 0 1 


64 = 0 0 0 


\ 1 0 0 



r 

1 

0 \ 

1 , 63 = 

1 

0 

0 

/ ' 


0 

V 

\ , 


0 

0 \ 

, 65 = 

0 

0 

1 

/ ' 


1 

0 / 


(37) 


Thus, any element in su(3)/so(3) can be written as a linear combination of 
these matrices, that is for any X G su(3)/so(3) there exist (cci,... ,x^) such 
that 

5 

X — ^ ^ 

n—l 

In particular, we say that any two elements X and Y G su(3)/so(3) are 
equivalent, and thus are elements of the same orbit, if there exists a matrix 
S{(pfi,tp}i,9fi) G 50(3) such that 


S{^r,iPr,9r) XnCn S \ipR, ^jjR, 9 r) = Y ynSn 

\n=l / n=l 

where x := {xi, ..., X 5 ) and y := (yi,..., j/ 5 ) correspond to X and Y, respec¬ 
tively. Similarly, this relationship can also be represented by a 5 x 5 matrix 
mapping x onto y. In particular, x ~ y if and only if there exist ipR, ij^R and 
9r such that 

Mc{‘fR,tljR,9R)x = y 

where the particular form of Me is given in Appendix ID. 41 Having estab¬ 
lished an equivalence relation that describes the orbits in case of the Cartan 
representation, we are now in a position to describe its orbit space. 

Lemma 20. Let X and Y be two real symmetric matrices. Then there exists 
a matrix S G 50(3) such that 5A5“^ = Y, if and only if the characteristic 
polynomials of X and Y agree, that is x(Al) = x(H). 


Proof. Omitted. □ 

Based on Lemma [SHI we deduce that the orbits described by the Cartan 
representation are separated by the coefficients of the characteristic polyno¬ 
mials of the corresponding matrix. In particular, the characteristic polyno¬ 
mial x{^) of S' matrix X G su(3)/so(3) written in terms of the basis E in 
dSil) is given by 

x{X){k) = — + k [xl+ X 1 X 2 +x\+x\+xl+ xl) 

— x\x2 — Xix\ + Xix\ — Xix\ -\- X2x\ — X2x\ -\- 2x^XiX^. 
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Thus, the two invariant polynomials for the group action written in terms of 
the Cartan representation are 

Ici '■=x\ + X1X2 + X2 + x'^ + x 1 + x'^ and 

Ic 2 ■= — x\x2 — X1X2 + X\x\ — Xixl + X2x\ — X2x\ + 2 x 3 XiX^. 

In particular, both of these polynomials satisfy the relationship 


Icj{Mc{(Pr,iPr,0r)x) = Icj{^) 

for j = 1, 2, respectively, and for all {^pr, tpR, Or) S [0, tt)^ x [0, 27r). Using the 
isomorphism given in (ED), the invariant polynomials that correspond to the 
group representation based on the Wigner rotation matrix in ID. 21 are given 

by 


T 2 I 2 I 2 I 2 I 2 j 

— X-^ “h X 2 X^ “h X^ X^ RrLQ 

^ _ 2x1x3 , ^_ , xlx3 ^2^ , 

12 —-7=—I" 2X1X2X4 H-- x^x^ + 


2 xi 


Xzxl 


73 


73 


373 73 73 


2X3X1 2 

“)“ X^X3 • 


In order to follow the programme outlined at the beginning of this section, 
it would now only be left to show that Ji and I 2 are generators of the ring 
of invariant polynomials. However, in this case, we already know that both 
polynomials separate the orbits by using the isomorphism and Lemma 1201 


4.2 Reduction to a two-dimensional orbit space 

Having found the fundamental invariants corresponding to the group action 
of SO(3) acting on we have an explicit description of all orbits as each 
orbit corresponds to exactly one constant value in the image of the two 
invariants. In order to reduce the state space to 5, we need to verify that 

{/i(a),/ 2 (a) I a e M'’} = {/i(a),/ 2 (a) | a G 5} . dSH) 

In particular, we can express the left hand side and the right hand side as 


00 


{/i (a),/ 2 (a) 

a G = [^ {{r, / 2 (a)} a G and /i(a) = r 


r—0 



00 

= U 7 } X 

min / 2 (a), max / 2 (a) 

/i(a)=r IilB.)=r 


r=0 

. aGK*^ aGR*' 

{/i(a),/ 2 (a) 

00 

1 a G 5} = M {{r, / 2 (a)} 1 a G 5 and /i(a) = r} 


r=0 



00 

= |J{r} X 

r—0 

min / 2 (a), max / 2 (a) 

Ilia)=r Ii{a)=r 

a£S a£S 
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respectively. Moreover, we observe that {/i (a) = r} is in fact a five-dimensional 
sphere which is connected and compact. Rewriting the problem like this, we 
observe that we are interested in the extrema of I 2 for points on the sphere 
with radius r. Since the sphere is compact, we know that these extrema 
are attained and and an application of the intermediate value theorem gives 
us that I 2 attains in fact all values in between its extrema. Taking these 
arguments into account, the verification of p6p has thus been reduced to 
checking that the extrema defining the intervals agree for all r > 0. Because 
both invariant polynomials Ii and I 2 are homogeneous, it actually suffices 
to verify this condition for the case r = 1. Moreover, since —/ 2 (a) = J 2 (—a) 
and a e if and only if — a S S^, it is enough to prove the equality of the 
maxima. Hence, we only need to show that the following relationship holds 


max 

aGK®,7i(a) = l 


max 


A direct calculation shows that in both cases 


max A (a) = max A (a) = —— 
aG 5 ^ ^ 3^3 


if /i(a) = 1, which thus proves the claim, and we conclude that the state 
space can in fact be reduced to the two-dimensional space S. 


5 Solving the reduced bifurcation equation 

Having shown that the equivariant state space of the bifurcation equation 
can in fact be restricted to the two-dimensional space S in (1351) . we are now 
in the position to solve this reduced problem. We denote the reduced form 
of the real bifurcation equation by which only depends on the two state 
variables oq and 02 . Its approximation up to fourth order is given by 


/reai_2(a>'^)(p) =0^ 

/reai_i(a,A)(p) =0, 



497r(l -b 327r)2 ^^5795328(1 -b 327r)2 



965888(1 -b 327r)2 



1931776(1 -b 327r)2 


/reaii(a,A)(p) =0, 
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/reai2(a,A)(p) =Xa2 + 


18 (Aa 2 aQ + Aa^) 
497 r(l + 327r)2 
9 5 


— 3207r)7r — l)a2ag 


3136(1 + 327r) 1792 

(567r(17 + 22407r) + 61)aiao 
241472(1 + 327r)2 


25872(1 + 327r)2 

/ 2 , 3n 

(0200 + 02 ) - -;^:77a2ao 


224 


(38) 


In particular, we observe that this equation can also be rewritten up to third 
order in Oi for i G {0, 2} and A as 


A) — (A + d(Oo + O2)) 


+ c 


Oq - “ 2 ^ 

— 2 ao 02 J 


(39) 


where c := and d := oTogiTfsSTry ~ r^j- However, we need to find 
conditions ensuring that we consider the equation up to a sufficiently high 
order, in order to guarantee that the characteristics of the bifurcation do 
not change. Such a problem is called recognition problem and it will be 
presented in the context of the two-dimensional bifurcation equation in (1391) 
in the following two sections. First we establish some theory in Section [5A] 
before this theory will be applied to our particular problem in Section 15.21 


5.1 The recognition problem in two dimensions 

The aim of this section is to derive non-degeneracy conditions which ensure 
a one-to-one correspondence between the solutions of the two-dimensional 
reduced bifurcation equation and an algebraic equation of a simple form of 
lower order. This problem is called recognition problem. In our presentation 
we follow [71 Chapter 9]. 

A one-to-one correspondence between the sets of solutions of two func¬ 
tions /(u, A) = 0 and (;(u. A) = 0 for u G M” and A G M exists if / and g are 
strongly equivalent. 

Definition 21. Two C°°-functions f and g are strongly equivalent, denoted 
by f -- g if there exist S' : M’" x K™ L(M’",K") and U : x ^ 

such that 


det S(u, A) > 0, 

dU 

(7(0,0) = 0, (u. A) is positive definite 

/(u,A) = S(u,A)g(C/(u,A),A). 

However, often it is not straightforward to verify strong equivalence for 
particular cases. Instead, it is easier if the problem inherits certain symmetry 
properties that can be taken into account. In case of the reduced bifurcation 
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equation in (l38ll . the problem is symmetric with respect to the group action 
of S3 acting on the space In particular, S3 acts on any element u G 

by 

giU — u' for all i G { 1 ,..., 6 }, 

where gi are defined as 

^3={5zGM2x2,f = l,..., 6 } 

_ f /l 0\ /l 0 \ 1 / -1 V 3 \ 1 [-1 

I \^o ’ \^o -ij ’ 2 V-v^ - 1 ^ ’ 2 - 1 ; ’ 

1 (-1 y 3 \ 1 / -1 1 

2 1 y ’ 2 1 J y 

Moreover, it is easy to show that the group action is generated by the elements 

52 = (J and 54 = ■_f) ■ (40) 

Based on this representation of the group action of S3 acting on M^, we 
direct our attention to the concept of S' 3 -invariance, S' 3 -equivariance and S3- 
equivariant equivalence. 

Definition 22. A mapping h : —?► M is S3-invariant if 

h{gu) = h{u) 

for all g G S3 and u G 

Definition 23. A mapping / : x M —>■ is Ss-equivariant if 

f{gu, A) = gf{u, A) for all g G S3 and u G 

Definition 24. Two S3-equivariant -mappings / : x K —>■ and 

5 : X M —> are equivariantly equivalent, denoted hy f ^ g, if there exist 
S' : R2 X R2 ^ L(R2,R2) and C7 : R^ x R^ ^ R^ such that 

det S(u, A) > 0, 

dU 

[7(0, 0) = 0, (u. A) is positive definite 

S(gu, A) = gS{u, X)g~^,U{gu, X) for all g G S3,u G R^,A G R 
/(u,A) = S(u,A)g(t/(u,A),A). 

The aim of Section [5^ will be to show that the reduced bifurcation equa¬ 
tion is in fact Ss-equivariant and equivariantly equivalent to its approxima¬ 
tion up to fourth order. In order to verify the equivariant equivalence, we use 
the following concepts and results. 

Theorem 25 [281 Theorem 1]. Any smooth S3-invariant mapping h can he 
represented as 

h{u) = h(hi(u),/i2(u)) 
where hi{u) := u\ + and /i2(u) := uf — Suiu^- 
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Theorem 26. Any S^-equivariant mapping / : x M —>■ can he repre¬ 

sented by 

/(u,A) =a(/ii(u),/i 2 (u),A) + 5 (/ii(u),/i 2 (u),A) 


with /ii(u) = Ui+ U 2 and /i 2 (u) = uf — Suiu^- 

Based on these results, we are now in the position to state the following 
proposition which yields the existence of a normal form under appropriate 
conditions. 


Proposition 27. Let / : x M —>■ be of the form given in Theorem 

Then 


/(u,A)^eA (“ 1 ) 


/uf - uiN 

\^-2uiM2 J 


for e := sgn(||) if and only if 


dcL 

a( 0 , 0 , 0 ) = 0 , &( 0 , 0 , 0 ) 7 ^ 0 and —( 0 , 0 , 0 ) ^ 0 . 


5.2 Existence of a transcritical bifurcation of the Onsager free-energy 
functional up to equivariance 

The aim of this section is to find non-degeneracy conditions that allow us to 
reduce the problem of solving the bifurcation equation in (15^ to an algebraic 
equation of a simple form and lower order. In particular, we show that /real 
is S' 3 -equivariantly equivalent to a mapping G : x K —>• by using the 
symmetry properties of /real- 

In order to prove that /real is ^a-equivariant, it is sufficient to show that 
it is equivariant with respect to the generators 52 and 34 of the group action, 
see go]). This in turn is a consequence of Proposition [191 Notice that for two 
particular choices of Euler angles we have the following two Wigner matrices 

/I 0 0 0 0 \ 

0 0 0 -1 0 
0 0 10 0 
0-100 0 
Vo 0 0 0 -1/ 

/ 0 0 0 -1 0 \ 

-1000 0 

0 0 -i 0 . 

0 10 0 0 
V 0 0 ^ 0 -i / 

Both linear maps leave the subspace S = {(0,0, cc, 0, j/)|a:, ?/ G R} invariant. 


Ri= R (f , 7 r, 7 r) 

and 

R-2 = 
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Merely viewed on the space S, they correspond to the two elements that 
generate S 3 . Hence we conclude that /real is S' 3 -equivariant and we know 
from Theorem [21] that it can therefore be written as 

/reai(u,A) = a(/ii(u),/12(u), A) -f 6 (/ii(u),(u), A) ■ 

Moreover, we can further deduce that it reduces to the normal form 

\U 2 j \-2U1U2J 
with e := sgn (ff) if and only if 

Od 

a(0, 0,0) = 0,6(0,0, 0) 7 ^ 0 and —(0,0, 0) ^ 0. 

oA 

Since it is not trivial to find the explicit form of the coefficients a(-, ■) and 
6 (-, •), it is not straightforward to verify that these conditions hold in case of 
the bifurcation equation in (l32ll . 

However, it is sufficient to show that they hold in case of the reduced 
bifurcation equation in (IHH because the first two derivatives agree. In 
this case the coefficients are given by 

S(/zi(a),A2(a),A) = 

Recall that Remark [HI shows that /real Is a fourth order approximation to 
/real- In /(^^i have dropped the fourth order terms so that it agrees with 
fs up to third order and therefore also verifies the non-degeneracy condition. 

Thus, we see that all three non-degeneracy conditions in Proposition 1271 
hold and we can conclude that there exists a one-to-one correspondence be¬ 
tween the set of solutions of in (|H 1 ) and 



The set of solutions of G can easily be computed. It is given by the trivial 


solution 




and the three branches 


^ ( 0^) 4 (ys) ^ (-^) ■ 

However, looking at the rotational symmetries of these three branches, we 
see that they in fact coincide after an application of a rotation. Thus, the 
solutions corresponding to all three branches are equivalent which proves 
the existence of a simple linear branch and thus the existence of a unique 
transcritical bifurcation up to rotation. 
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6 Uniaxial Solutions 

Recall that the Euler-Lagrange operator corresponding to the Onsager tree- 
energy functional is given by 


E{cj^,X) -xm- 


1 




k{p ■ q)e-Kp{-(j){q)) dq. 


© 


'§2 


Now, we restrict our attention to the uniaxial solutions of this problem, that 
means that we assume that any solution ps : L^(S^) is axially symmetric 
with respect to the z-axis. Thus, writing p in terms of spherical coordinates 
9 G [ 0 , 7 r) and tp G [0,27r) (see Appendix [Xl for details), we assume that it 
only depends on 9 and is independent of ip. In this case, where 


dp = 


dp 


sin 6 * d9. 


/S 2 


the Onsager free-energy functional in m restricted to the set of axially sym¬ 
metric probability densities can be rewritten as 

/>7T 

Es{ps) ■=kBT I ps{9p)ln{ps{9p))sin9p d9p 


A 2 


K{9p, Pp, 9q, Pq)ps{9q)ps{9p) sin 6 »p siu d9pd9qdppdpq 


for A := [ 0 , 7 r] x [0,27r]. The Euler-Lagrange equation corresponding to this 
restricted problem does not differ from the one of the full problem. The 
reason is that the Euler-Lagrange equation is derived by computing 


-j-J' s(ps + X) 
de 

= [ Zs(9p) 


e=0 


X In Ps (9p) -|- / AT(9pj Pp^ 9q<j Pq^Ps {9q^ sin 9q d9qdpq 


d9„ 


where Zg ■ [ 0 , 7 r) — >■ M such that Zs{9p) sin 9p d9p = 0. In particular, 
the term within the square brackets does not depend on the variable pp, for 
details see LemmaH^in Appendix!^ This allows us to apply the fundamental 
lemma of the calculus of variations with respect to the variable 9p and gives 
us the Euler-Lagrange equation 

Alnps(6>p)-f- / K{9p,Pp,9q,Pq)ps{9q)sm9q d9qdpq = -XinZs 

Ja 


where A := ksT and 




exp 1-1 


r2TT 


K{9p,pp,9q,pq)ps{9q)sin9q d9qdpq sin9p d9pdp^^ 


as before. Introducing the thermodynamic potential (ps ■ 

7 

IA 


as 


A 


<(>77) : 


K{9p,Pp,9q,Pq)ps{9q) Sin 9q d9qdpq , 
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it follows that Ps{dp) = ^ exp(—(/)s(0p)) and we can rewrite the Euler- 

Lagrange equation as 

Ms{0p)-^ I K{0p,(pp,0q,(pq)exp{-(j)s{9q)) deqdipq=0. 

Thus, we conclude that the Euler-Lagrange operator for uniaxial probability 
distributions is equivalently given by 

Es{(j),X)=X(j)s{0p)-^ K{0p,tpp,0q,(pq)exp{-(l>s{0q)) d0qdipq. (41) 

A well-known concept that establishes a relationship between the symmetric 
solutions of a minimisation problem and the symmetrised problem is the 
principle of symmetric criticality m- It states that each critical point of the 
symmetric problem is a symmetric critical point of the general problem. In 
the above setting, this principle can be verified explicitly. Since the Euler- 
Lagrange operators Eg in dm) and E in d3 have exactly the same form, we 
can consider the corresponding Lyapunov-Schmidt decompositions for the 
general and the symmetric problem simultaneously. In particular, we split 
(j) = u + V and /ps = Ug Vg, respectively. The first step consists in solving 

(1 — P)E{u -I- z;. A) = 0 and 
(1 - P)Eg{ug + Vg, A) = 0 

for V in terms of u and Vg in terms of Ug, respectively. Taking the derivative 
with respect to v and Vg, respectively, we see that we can use the implicit 
function theorem in order to solve uniquely for Vg and v. Because E(ug + 
Vs,X) = Eg(ug -|- zzs) A), it follows that 

v{Ug) = Vg{Ug) 

is axially symmetric. Therefore, solving the bifurcation equation of the full 
problem with a _2 = a_i = oi =02 =0 yields all uniaxial solutions. In 
particular, this equation is given by 

ISAag ^ yf(4487r(13 + 28007r) + 577)a4 

’ 497 r(l-b 327r)2 5795328(1 -b 327r)2 

f 9 5 A 3 1 2 

V 3136(1-b 327r) 1792 J 448 

This is now a bifurcation problem with a one-dimensional state variable 
and we say that it undergoes a transcritical bifurcation at ( 0 , 0 ) if ( 0 , 0 ) is a 
non-hyperbolic fixed point, that is 

df 

/( 0 , 0 )= 0 and^( 0 , 0 ) = 0 , 
oaQ 


and if the non-degeneracy conditions 
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hold (these conditions can be found in the literature, see [37l (3.1.65)-(3.1.68)]) 
It is easy to see that all of these conditions hold for /s(ao,A) and we may 
therefore deduce that a transcritical bifurcation occurs. Since we do have 
a unique transcritical solution in case of the full problem and the symmet¬ 
ric problem, we conclude that these are the same. Therefore we proved the 
following statement. 

Theorem 28. All solutions of the Onsager free-energy functional are uni¬ 
axial in a neighbourhood of the trivial solution fi = 0 and locally around 

A = A2. 

7 Qualitative behaviour of the global bifurcation diagram 

Finally, we complete our bifurcation analysis of the Onsager model by proving 
some properties of the global bifurcation diagram. In Section lTTl we show that 
any minimiser of the free-energy functional in ([ 2 ]) is bounded and that we can 
therefore restrict our attention to the set of bounded probability densities. 
Using this result, we show that the free-energy functional is in fact strictly 
convex for high temperatures and that its trivial solution, po = -^ is the 
unique global solution. In Section [7j2] we prove that the trivial solution is a 
local minimiser for high temperatures and that it is not a local minimiser 
for low temperatures. Finally, we show in Section 17.41 that any bifurcation 
branch either meets infinity or that it meets another bifurcation branch, thus 
proving that all bifurcation branches are continuous and do not end suddenly. 

7.1 An upper bound on all admissible probability density functions 

By constructing a test function p* that is bounded above by a constant C*, 
we prove that iF{p,X) > iF{p*,X) for all admissible probability densities p 
and thus we show that we can restrict our attention to probability densities 
that are bounded above. This result holds for all interaction kernels that are 
continuous and symmetric, see Assumption [5] (a)-(b). 

Lemma 29. Let p G L^(S^) be a probability density and suppose that there 
exists a set A C of positive measure such that p{p) > C* := exp(16M/A) 
for all p G A where M := max^ ^^§2 K{p,q). Then there exists a modification 
of p denoted by p*, such that p* is bounded from above by C* and J^{p*, A) < 


J-(p,A). 

Proof. Let 

p*{p) := C'*l{p|p(p)>c*} + P(7')l{p|C*>p(p)>^} + (pip) + T')^{p|p(p)< 


where 7 := 
following notation 


/s2(p(p)-C*)l{p(p)>c«} dp 

I{P|2^>P(P)}I 


. For the ease of presentation, we use the 
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C := Ip 


2tt 


> pip) 


It is easy to see that IC*! > 0 since otherwise pip) dp > 1. By the choice of 
7 , we make sure that p* is a probability density that still integrates to one. 
Observe that C* = exp(16M/A) > 1 > ^ and hence A 0 B = 0. In order to 
prove that 


F{p,\)-Fip\\)>Q, 

we deal with the entropic and the interaction term separately. We begin by 
deriving a lower bound on the entropy term. In particular, 

Ti{p,\)-W,\) 

=A [ p{p) \n{p{p)) dp- X [ p*{p) ln(p*(p)) dp 

=>' [ [fipip))-fiC*)]dp +X [ [fipip))-fipip)+7)] dp 
J A JC 

where fix) := a;ln(a;). Moreover, defining g(p) := p(p) — C*, we may apply 
the fundamental theorem of calculus (notice that the following claim holds 
trivially for ^(p) = 0, so we may assume that gip) 0 for all p G S^). Hence 

J-i(p,A)-J-i(p*,A) 

f /to(p)+C-)-/(C») f /Wp)+7)-/W) ^ 

Ja dip) Jc 7 

=A f fifi)gip) dp- Xj [ /'(C2) dp 
Ja Jc 

where ^i(p) G [C*,gip) + C*] and f 2 ip) G [p(p),p(p) + 7 ] with p € C. Since 
fix) = ln(x) -I- 1, which is monotonically increasing, it follows that 

.Fi(p,A)-J-i(p*,A) 

>A(lnC'* + I) / gip) dp-jX [ (ln(p(p) -h 7 ) -b 1) dp 
Ja Jc 

>A(lnC'*-bI) J gip)dp- 7 X J ^In 

where the last line follows from the fact that p G C and hence p(p) < ^ for 
all p G C. 

We now claim that 7 < ^ • Since p is a probability density, it follows that 
/^(p(p) — C*) dp < 1. Hence 

Is^ipjp) - Cfl{p(p)>c*} dp ^ 1 

\{p\^>Pip)}\ ~ \{P\^ > Pip)}\ 

and we only have to show that \C\ = |{p|^ > pip)}\ > 27r. Suppose the 
contrary, that is \C\ < 2tt and thus |H U H| = |{p(p) > ^}| > 2tt since the 




-I- 1 dp 
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total measure of the unit sphere is given by 47 r and thus \A\J B U C\ = 47 r. 
Moreover, since p is positive, 


/ p{p) dp> pip) dp > 

Jqc 


1 

27r 


'C‘ 


1 dp > — 27r = 1 
27r 


contradicting the fact that p integrates to one. Therefore 7 < 


Observing that gip) > 0 for all p G A and using the monotonicity of the 
logarithm as well as the facts that 7 < In ^ < —1 and A, 7 , \C\ > 0, we 
deduce that 


J-i(p,A)-J-i(p*,A) 

>A(ln C* + 1) f gip) dp — X 

J A 

>Aln(C'*) [ gip) dp. 

J A 

This gives us a lower bound on the entropic term and we can direct our 
attention to a lower bound on the interaction term. Using the definition of 
p*, it follows that 


ln|^ + 5^l+l)7|C| 


7^(aA)-7^2(p*,A) 

“2 
_1 
“2 


Kip. q)[pip)pi<l) - p*ip)p*iq)\ dq dp 


/S 2 xS 2 


' Kip,q)hip,q) dq dp-— 
AxA ^ 


■ 7 


Kip,q)pip) dq dp - 


/ Kip,q)dqdp 
ICxC 

Kip,q)piq)gip) dpdq 


ICxC 


f AxB 


Kip, q)piq)gip) dq dp - -fC* 


Kip,q)piq) dq dp 


JAxC 

7 / Kip,q)piq) dq dp 
■JbxC 


lAxC 


where we used the symmetry of the kernel and again we defined gip) := 
pip) — C* and hip,q) := pip)piq) — C*C* which are both non-negative for 
p,q € A. Assuming without loss of generality that the interaction kernel is 
positive, we observe that all terms with a positive sign are positive while 
only those with a negative sign are negative. Neglecting the positive terms, 
we may therefore deduce that 


K2ip,X)-Mp\^) 


>- — Kip,q) dqdp--f 

^ JCxC 


Kip, q)pip) dq dp 


CxC 


-7 


Kip,q)pip)dp] dq-"i 


Kip,q)piq) dqdp (42) 
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> — (47r)^ — — AnjM — AttM^ 

> — 7^M87r^ — 12tt^M 

where we used the fact that p{p) > C* for all p S ^ in (11^ . Since 
J-2{Pj A) — J-2{p*j A) > — Att^M — VItt^M = — 167r7M, 
we obtain a lower bound on the interaction term. 


Thus, we are now in the position to derive an overall bound on T. Putting 
both lower bounds together and using the definition of 7 = 
obtain that 


A) — A) >A ln(C'*) j g{p) dp — 

J A 

=Aln(C'*)|C'| 7 - 167 r 7 M. 

Finally, making use of the assumption that C* = exp(16M/A) and using the 
fact that \C\ < 2tt, the statement follows because 

J-(p,A) - J-(p*,A) >Ai^|C|7-167r7M = 16M7(|C'|-7r) >0. 

A 


□ 


Corollary 30. Let two sets A and B be given by 

.A := —7> R| J p{p) dp = 1, p > 0 

S := |p : —>■ M| J p(p) dp = 1 , p > 0 , p < exp(16M/A) 

Then any solution of the minimisation problem min^J^(p, A) is in fact a 
solution of the redueed minimisation problem ming J'{p, A). 

Proof. The statement follows directly from Lemma [29l □ 


7.2 Local properties of the trivial solution po = ^ 

We concentrate on proving conditions which ensure when the trivial solution 
is a local minimum and when it is not. 

Proposition 31. For all A > A 2 , the trivial solution po = ^ is o local 
minimum. For all 0 < X < X 2 , the trivial solution po = -^ is not a local 
minimum. 
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Proof. The second variation of the Onsager free-energy functional in ([S]) in 
L°° is given by 


52 


I{z,X) = DlTip) = i:^J^{p + ez 

= A 




/S 2 


pip) 


de^' 
dp - 


/S2xS2 


e=0 

k(p ■ q)z{p)z{q) dpdq 


(43) 


with 2 ; G L°°(§2) such that Jg z{p) dp = 0. The Stone-WeierstraB theorem 
guarantees that we can find Om.i such that 


J2am,iYr-z 


< e 


for some e > 0. Using the fact that the interaction operator U turns into a 
multiplication operator when it is applied to a spherical harmonic, 


u{Yn = 

we can rewrite the second variation as 


Pi if I is even 
0 if Hs odd, 


00 I 


00 I 


I{z,X)=AttX / (y] ^ aimYf^ip)] dp 

n / 00 Ip \ L 

^ LAY T. aim lk{p-q)Yf^{p)dp ]YY aimYriq) dq 

p / 00 I \ ^ 

=47rA ^^irnYrip) 1 dp 


1—1 m— — l 


Y Y aimPiYriq) YY aimYriq) dq 


I i 


IS even m 




1—1 rn——l 


where we skip the term Z = 0 so that p + z still integrates to 1. Viewing the 
spherical harmonics as an orthonormal basis in r{S'^), it follows that 

00 Z 00 z 

I{z, A) = 47 rA y] y \aim\^+ y y pi\aim\Y 

1—1 m— — l I is evenm= —Z 


Hence the second variation is positive, that is I{z,X) > 0 for all 0 7 ^ z £ 
r{S^), if and only if 47 rA > —p 2 - Thus, the trivial solution po is a local 
minimiser of the bifurcation equation for all 
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Moreover, we observe that the following holds 

=-^ (s)+ (3;) + k" a I -(<•)+ 


^ /(^(p)^-^), + h-o.t 

<0 if A<A2 


= j' 


^ j + 

If we choose z close enough to zero, then the higher order terms are negligible 
and we see that p = ^ is not a local minimiser for all A < A 2 . □ 

7.3 Existence of a unique solution for high temperatures 

We prove that the Onsager free-energy functional is strictly convex in p and 
therefore that the isotropic state, represented by the uniform distribution 
Po = is the unique solution for high temperatures. 

Theorem 32. Let k(-) be a continuous interaction potential and let A be 
given such that the following inequality holds 

SttM exp(16M/A) < A. 

Then the Onsager free-energy functional T is strictly convex on the set of all 
probability densities p bounded by C*, see Lemma 

Proof. Recall that the second variation of T is given by 


A 


\v) 


dp - 


k(p ■ q)z(p)z{q) dp dq > 0, 


JS'^ P{P) JS2xS2 

see (H5)) . We know from Lemma HH] that it is enough to consider the min¬ 
imisation problem among the set of probability densities p € L^(S^) which 
are bounded by a constant C* = exp(16M/A) where M denotes an upper 
bound on the kernel k{p ■ q). Therefore an application of Jensen’s inequality 
and imposing the condition that 8 ttM < X/C* yield 

zHp) 


A 


/S 2 


P{P) 


dp - 


/S^xS^ 


k{p ■ q)z{p)z{q) dp dq 


Js^ 

zHp) 

pip) 

dp — 

M 

[ \zip)\ 

Js^ 

zkp) 

pip) 

dp — 

MAtt / 

z{p)^ dp 

JS2 

z\p) 

pip) 

dp — 

— [ 

z{p)‘^dp 

1 

< 0 

1 

Al 

/ 

ip) dp 


> 0 , 






dp 

Jtt 
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proving that the Onsager functional is strictly convex. Rearranging this con¬ 
dition and plugging in the constant C* = exp(16M/A) yields the above 


inequality SttM exp(16M/A) < A. 


□ 


Theorem ED proves that the Onsager free-energy functional is a strongly 
convex functional for all probability densities p G S whenever the temper¬ 
ature parameter satisfies the inequality stated in Theorem Numerically, 
this value can be approximated by A* ~ 38.205 when M = 1 (for example in 
case of the Onsager interaction potential). 

Corollary 33. Let fc(-) be a continuous interaction potential and let A* be 
chosen such that SttM exp(16M/A) < A. Then the uniform density po = ^ 
is the unique solution to the Onsager free-energy functional for all A > A*. 

Proof. The uniform density po = is indeed a solution to the Euler- 
Lagrange equation for all A S M. By Corollary EOl and Theorem EH the 
necessary optimality condition is also sufficient PI Theorem 5.2.4]. □ 

7.4 Continuity of all bifurcation branches 

We conclude this section by showing that every bifurcation branch of the 
Onsager functional is either unbounded and meets infinity or that it meets 
another bifurcation branch. 

The Euler-Lagrange equation in terms of the thermodynamic potential is 
given by 



which is of the form Aid — U where U is compact, see Lemma PI Therefore 
we can rewrite the bifurcation problem as 


lc/(</.,A) = </> 


and it is thus of the form of bifurcation problems considered by Rabinowitz 


[13- 


Theorem 34 pi Theorem 1.3]. Let G : E x ]S. ^ E be a compact and 
continuous non-linear operator and consider the problem 

u = G{u, A) 

where u G E and A £ M. Moreover, let 

G{u, A) = \Lu -\- H{u, A) 

where Ll{u,X) is 0(||u||) for u near the origin uniformly for bounded A and 
let L he a compact and linear map on E. If X is a real non-zero eigenvalue of 
L of odd multiplicity, then the solution set {(u. A) \ u = G{u, A)} possesses a 
maximal connected and closed subset Gu such that (0, A) G Cu and Gu either 

1. meets infinity in E, or 
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2. meets (0, A) where X ^ X is a real non-zero eigenvalue of L. 

In particular, Rabinowitz shows that any eigenvalue of odd multiplicity is 
a bifurcation point and that any branch bifurcating from the trivial solution 
must either meet infinity or meet another bifurcation branch. In case of the 
Onsager functional, all eigenvalues have multiplicity 2n -|- I, thus we know 
that (0,As) is not an isolated solution but that it is instead a member of a 
non-trivial closed connected set described by Theorem 1341 


8 Conclusion 


We found a general method to derive the eigenvalues of integral operators of 
the form 



K{p,q)p{p)piq) dpdq 


depending on kernels satisfying Assumption [S] and we derived an explicit 
expression for the eigenvalues of the Onsager free-energy functional involv¬ 
ing the Onsager kernel. Based on this result, we know from that all 
eigenvalues are bifurcation points and we proved that a transcritical bifur¬ 
cation occurs at the temperature r* = X 2 /kB where A 2 = ^ denotes the 
largest of these eigenvalues. Moreover, we showed that the corresponding so¬ 
lution is uniaxial. Globally, we showed that the trivial solution p = ;^ is the 
unique solution for high temperatures. Finally, we verified that all bifurcation 
branches either meet infinity or they meet another bifurcation branch. 

In order to obtain a complete bifurcation diagram, it is left to show that 
no isolas occur, that means that we do not have any bifurcations that are 
isolated from all other solutions to our problem at hand. Our approach is not 
at all limited to the Onsager kernel and not even to quadratic forms given 
by integral kernels. This will be demonstrated in forthcoming work. 
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A A brief introduction to spherical harmonics and associated 
Legendre polynomials 


The spherical harmonics are defined as 

Yr{^,e) -Nime'-^^pr{cos 6) 

where ip € [0, 27r) denotes the polar angle and 9 € [0, tt) the azimuthal angle. 
The subscript I denotes its degree while denotes its order. Nim is a 

normalisation constant given by 


Nim := (- 1 )" 


2Z-Piy/^ yz-m)!y/^ 


Itt 


(Z-bi 


and PP is the set of associated Legendre polynomials given by 

|?T=rC-i)‘ 


One of the most important facts of this set of spherical harmonics is that they 
build an orthonormal basis of L^(S^). In general, spherical harmonics are known 
as eigenfunctions of the Laplacian. However, it seems that they can be understood 
in a much broader sense as eigenfunctions of rotationally symmetric operators. 

For m ore infor mat ion on spherical harmonics and associated Legendre polyno¬ 
mials see [33] and |34 |. 


B Technical results used in the derivation of an explicit 
expression for the eigenvalues of the Onsager functional 


Lemma 35 (Supplement to Theorem jS]). 


r 


X 


E 


(2Z-bl)r! 

2(’-0/2(l(r-Z))!(Z-br-bl)!! 


Pi{x). 


(44) 


Proof. The claim follows from an induction based on the reccurrence relation for 
Legendre polynomials 


Z -t- 1 
2Z -b 1 


Pi+i{x) -b 


Z 

2Z + 1 


xPi (x) 


Pi-i{x) 
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with the convention that P_i = 0 [33]. The base case r = 0 is trivial because 
Pq = 1, so let us move on to the actual induction step. Multiplying both sides of 
(Ifill with X yields 




E 

l = r,r — 2,.. 


{2l + l)r\ fl + 1 

2(r-0/2(i (r - + r + 1)!! V2^ + 1 


Pi+i{x) + 


I 


21 + 1 


Pi-i{x) 


Changing the index in the two parts of the sum corresponding to Pi+i and Pi-i, 
we obtain 


X 


r+1 


= E 


/r!(/+ r + 2) + 2i(r + 1-/)(/ +l)r! 
2C+i-0/2(i(r + 1 - l))\{l + r + 2)!! 


+ 


r + 1 


(; + l + r)!! 1 


++1 


(x) + 


_3d_ 1 Pn(r) 

_5H_ 2 p / \ 

2('--2)/2(i(r-2))!(r+3)!! 


if r is odd 
if r even 


Notice that we excluded the boundary cases from the overall sum. Putting every¬ 
thing together yields 


X 


r + 1 


= E 

i = r + l,r— 


(2Z + l)(r + l)! 

2C+i-')/2(i(r + 1 - l)y.{l + r + 2)!! 


Pi{x) 


which concludes the induction step. 


□ 


Lemma 36 (Supplement to Corollary IIUII . The change of integral and infinite 
sum in (I13|) i s valid in case of the Onsager kernel. In other words, we can verify 
condition (1161) . stating that 


^^(4Z + l)^/+r|<oo, 

1=0 r=l 


for 


( _ 4n(2r)\^ _ if I < r 

(Z ■— .1 (l-2r)(r!)2(4'-)2>--'(7--i)!(2i+27--)-l)n V ^ ' 

■ \0 ifl>r' 

Moreover, the eigenvalues of the interaction operator equipped with the Onsager 
kernel decrease faster than l~^. 


Proof. With regard to m, the Onsager kernel written in terms of spherical har¬ 
monics is given by 


koip-q) = J2J2 E 


47r(2r)!2 


+ U .(1 - 2r){r\y{M)2r-i{r - l)\{2l + 2r + 1)!!^ 


Y27{p)Y^2i{q). 


da 
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The coefficients c [ ar e all negative except for Cq = dvr. Taking this into account and 
applying Lemma [37l we know that 


1=0 r=l 

oo oo oo 


r=l 

c 

',3/2 


——^(4z + -I-1 — cs -b 2c! 


+ ra + i)ra + 2) +l 2r(i)r(2) 


= 9Z.(4' + i) r(; + i)r(/ + 2 ) + I”" +®"j- 


(45) 


Using Gautschi’s inequality m, 

r{x + i) 


r{x + a) 

it follows that 

r{i + ^)r{i-i/2) ^ rii-1/2) 


< (a: -b 1) “ for any 0 < a < 1, 


< 


r(z-i/ 2 ) 


r(z-b i)r(z-b 2) (;-bi)i/2r(z-b2) {i + i)^/^rii)iii + i) 

1 ^ 1 
;(; + i )2 ^ ¥' 

Inserting this result into (|45l) yields 

', 3 / 2 , 7 -| , TT (4^ + 1 )^/^ 


(46) 


OO OO 


y:y:(4 z+i)^/^icn<^E' 


Z3 


1 = 0 r=l 1 = 1 

Comparing this with the harmonic p-series, we observe that 

1 


-b . 


^ - (4Z + 1)3/2 

2 ^ 


1=1 




+ + 8^) ~ C'l p 72 +C 2 <00 


for some constants Ci and C 2 which therefore proves the first part of the claim. 
The second part follows easily from (l46l) . □ 


Lemma 37 (Supplement to Corollary I lOL Let s £ 2N. Then 


OO 


E 

r=sl 2 


47r(2r)!^ 

(1 - 2r)(r!)2(4^)2’-V2(r - 


s/2)!(s-b2r + 1)!! 


^r(s/2+l)r(s/2-l/2) 

2r(s/2-b l)r(s/2-b2) 
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Proof. First of all, let us rewrite the sum by translating it by s/2. This yields 
47r(2r)!2 

2r)(r!)2(4'-)2’'-«/2(r - s/2)!(s + 2r + 1)!! 

_ 7r((2r + s)!)^ 

^ ^ (1 - 2r- - s)((r + s/2)!)223’-+'>-2r!(2s + 2r + 1)!!' 

In case of odd integers n = 2A: — 1 for A: > 1, the double factorial n\\ can be replaced 
by the expression n!! = | 35[ page 823]. Hence we can rewrite the terms such 

that 

7r((2r + s)!)^ 

^ (1 - 2r - s)((r + s/2)!)2(23’-+“-2)r!(2(r + s + 1) - 1)!! 

^ 7r((2r + s)!)"2-"’~+"(r + s + l)! 

(1 — 2r — s)((r + s/2)!)2r!(2r + s + 1)! 

_ / (2r+ s)! y (r + s + 1)! (2r + s - 2)! 2"^’'+^ 

\ + • 5 / 2 )!/ (2(r + s + 1))! (2r + s — 1)! r! 

Notice that we also expressed the first factor in the denominator as fraction of two 
factorials. This s um can now further be manipulated by rewriting it in terms of 
gamma functions |35[ page 1137], 


r=s/2 ' 


r(s +1) = s!. 


It follows that 


^ / (2r + s)! y (r + s + 1 )! (2r + s - 2 )! 2 ~^^+^ 

V (’'+ s/ 2 )!/ ( 2 (r + s + 1 ))! ( 2 r + s — 1 )! r! 

^/r(2 r + s + l)y r(r + s + 2) P{ 2 r + s - 1 ) 

/"(»■ +s/ 2 + 1)7 r( 2 (r + s) + 3 )) r( 2 r + s) r! 

^ / (2r + s)r(2r + s) ^ (r + s + l)r(r + s + l)r(2r + s - 1)2"^"+^ 
V (’' + s/2)r{r + s/2) J (2(r + s + l))r(2(r + s + l))r(2r + s)r! 

_ / r( 2 r + s) y r(r + S + 1 ) T( 2 r + s - 1 ) 2 "^"+^ 

^ ^f.\r{r + s/ 2 ) J r(2(r + s + l)) r(2r + s) r! 


Applying the duplication formula for the gamma function |35l page 1139] 

r(2z) _ 2^^-^r{z + ^) 

Ti^- ^ ’ 


to the first two factors and its rearranged form 
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to the two terms forming the third factor, we deduce that 

-^/r(2r-bs)y r(r-bs-bl) r(2r-b s - 1) 2 
+ r(2(r + s + l)) r(2r + s) 

'222^r(r-b s/2 + i) V 2i-2(’'+«+i)0F 

r=0 \ 


-2r+4 


= -E 




r(r -b s -b 3/2) 

22r+s-2p(^ -b sl2)r(r -b s/2 - 1/2) 


= -v^E 


22 ’-+=-ir(r-bs/2-bl/2)r(r-bs/2) r! 
r{r -b s/2 -b \)r(r -b s/2 - 1/2) 1 


rl' 


r(r + s + 3/2) 

Let {x)r denote the Pochhammer symbols which are given by 

(a;)o = 1 

-r( 2 ; -b r) 


(®)r = 


r{a 


for any x £ C and r G E. Rewriting our expression in terms of Pochhammer symbols 
[3], we obtain 


-v^E 


r(r -b s/2 + |)r(r -b s/2 - 1/2) 1 
-r(r -b s-b 3/2) r! 


V¥r(s/2 -b |)r(s/2 - 1/2) ^ (s/2 -b |)r(s/2 - 1/2), 1 


r(s-b3/2) 


E- 


(s + 3/2)r 


This sum can now easily be evaluated using generalised hypergeometric functions 

Bi- 

Definition 38. The generalised hypergeometric function pFq is defined as 

17 / N (a;i),(a;2)r ... (xp), z'’ 

pFq{xi,X 2 , ■.. ,Xp-,yi,y 2 ,... ,yq-, z) = ^ , . , . - 

(j/l)''(y2)r ■ • • (i/g)r T. 

where (x), denote the Pochhammer symbols. 


In particular, we are now in the position to use Gauss’s Hypergeometric Theo¬ 
rem which states that 

Theorem 39 (Gauss’s Hypergeometric Theorem, |35L page 1438]) 


2Ti(a,6, c; 1) = 


r{c)r{c — a — b) 
F{c - a)r{c - b) 


for TZ{c — a — b) >0. 

An application of Gauss’s Hypergeometric Theorem yields 


V¥r(s + |)r(s-i/2) 
r(2s-b3/2) 


oo 


E 


(s -b ^)r(s — 1/2), 1 
(2s -b 3/2), H 


V¥r(s+^)r(s-l/2) 

r(2s-b3/2) 


2 Ti(s + l/2,s-l/2;2s-b3/2, 1) 


V^r(s -b i)r(s -1/2) r(2s -b 3/2)r(3/2) 
r(2s-b3/2) r(s-b i)r(s-b 2) 

^r(s+|)r(s-i/2)r(3/2) 

F{s -b l)-r’(s -b 2) 
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Using the fact that r{ 3 / 2 ) = the claim follows. 


□ 


C Technical steps in the derivation of the bifurcation equation 

C.l Regularity of critical points and differentiability of E 

First we show that the interaction operator is a compact operator on the space 
for all s G R. 

Lemma 40. The interaction operator U is a compact operator —>■ 

for all s € K. 

Proof. By (I17II we see that U acts as a multiplication operator applied to the 
orthogonal basis of spherical harmonics. Due to Assumption [S] (d), we know that 
its eigenvalues converge to zero. Thus, 17 is a compact operator. □ 

The following Proposition shows that solutions to the Euler-Lagrange equation 
are smooth. 

Proposition 41. Suppose that (f G L°°(S^) satisfies 

= (Xs + X)(f>{p) - [ ko{p ■ dq = 0 . 

Z(cp) ys2 

Then </> G 

Proof. We rearrange the Euler-Lagrange equation as 

{Xs + X)cf){p) = f ko{p ■ dq. 

Z(<p) Js^ 

Notice that 0 is bounded, see Lemma M Thus, taking the limit for a sequence 
Pn —^ p on the right hand side, we see that it is continuous due to the dominated 
convergence theorem and so is the left hand side. Because ko{p- q) is differentiable 
almost everywhere, except for \p ■ q\ = 1, with an bounded derivative, we can 
swap integration and differentiation and conclude that 0 is differentiable. 

Having proved that tj) is continuous and differentiable, we want to show that 
also the derivative of (j) is continuous and differentiable. This will then allow us to 
proceed by an induction argument. 

By dU one way to define the C^-norm is 

ll<?!>llci := II^11 CO + , max 

l<t<j <3 


where 



Xi cost — Xj sint, Xi+i, ..., Xj-i,Xi sint -|- Xj cost, Xj+i,..., Xnj 

are the intrinsic derivatives on 5'"'“^. Here we consider the special case n = 3. An 
elementary calculation shows that 


x,jYr(p)vr(g) = -Yr(p)x.,Yr(q). 
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Thus, using the expansion of Ko{p,q) = see (IT8|, 

we can calculate 

[ Ko{p ■ dq 

Z[<p) Js^ 

^ p oo I 

oo I 

= EE^* 

1 = 0 rn= — l 
oo I p 

= E E / (-XijYr{q)e-^^‘‘^) dq 

in — / \ / 

1 = 0 rn= — l 

oo I p 

= E E / Yr{q)(-Xij<l,{q)e-^^‘‘^) dq 

in - / V / 

1 = 0 m= — l 

= ko{p ■ q) dq. 

In particular (As -|- X)Xij(j>{p) is continuous and differentiable. By iteration of this 
argument, the statement is proved. □ 


/ Xi^Yr{p)Yr{q)e-^^‘‘^ dq 
Js^ 


The Sobolev spaces on the sphere can be defined in terms of several equivalent 
ways through the coefficient in the expansion in spherical harmonics or in local 
coordinates p and 9. In terms of / = o-imYi^ the Sobolev norm can 

be defined as 

OO I 

Il/llrf^(s2) :=E E (47) 

1=0 m= — l 

In order to write down the equivalent norm in local coordinates, we follow [16] 
and use spherical coordinates p and 9 on the sphere, see Appendix]^ Hence the 
Riemannian metric in terms of local coordinates is given by 




( sin^ 9 0 

I 0 1 


and the Christoffel symbols are given by 

rk = 0 ,rk = rk = cote ,rk ^Q,rk = -\ sin 29,rk = rk = rk = 0. 

The 77® (S^) norm is then given by [!§] 

Il/llir»(s2) = E (/ Mg)) 

i=0 '' 

where denotes the covariant derivative. The corresponding norm is given by 

Subsequently, we are only interested in (S^) and thus we only need V/ and 
V^/. These are given by 

Xf = df 

{x^f).. = dijf-rkdku. 
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Hence explicitly, the norm can be expressed as 

||/|lH^ 2 ('g 2 ) = f sin edtpde^ + (sm~'^ e {d^ff + {deff'j sinOd^pde^ 

+ (/ sin”'^ 9 (^d<p<pf - i sin26»(9e/j + {d<p<pff (48) 

+2sin“^ 9 (dipipf) — i sin 29dgf^ sm9difid9'^ 

By the Sobolev embedding theorem m, 

II/IICO <C||/||^2(S2). (49) 

Lemma 42. E is infinitely many times Frechet differentiable as a mapping from 
(S^) (S^) . 

Proof. In order to check Frechet differentiability of E we need to verify that 

\\E{<p + ri) — E {<p) + Dr,E{(f>)\\^2(s'^\ 

lim -j—r- = 0. 

II’)IIh2(s2)“^° II’7|Ih2(s2) 

First we compute the first and second directional derivatives of E 

E (<(>) (p) = (A2 + A) <f>{p) - ^ J k{p- 

deE (0 + £77) (p) = (Aa + A) p(p) - ^ -—- 77^^ [ Hp ' q)e~’^^'^Em(i) 

dlE {(j> + eri) (p) 




Ti 


f e-^^‘’^-^’’<~‘^'>p^{q)dq (/ e-‘^<-‘>^-^^<-‘>^dq'^ ^ 

j k{p ■ q)e 

' (J e-'#>('i)-'’)('J)d(7)‘‘ J 

T2 

(50) 

(^f e-^‘-'^^-‘^<-'^'>piq)dqy 2 f dq 

y k(j) ■ q)e~'^^'^'’~'^'^iFi 






Ti 
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Thus, by the integral form of the remainder in Taylor’s theorem in e we have that 
E{4> + ri)-{E {<t>) + D^Em = + er,) (1 - e)de. 

The ||•||^ 2 ^s 2 ^-norm of the left hand side can be bounded by 


/ 

JO 


and our aim is now to establish an overall bound of the form 

||i? {4> + v) — {E {<f>) + DnE{(l>))\\f^ 2 (g 2 'j < ^(Il'?^ll/j- 2 ('s 2 ) , ||»7||^2('s2p | 




where 6 : R is locally bounded. This would yield that E is Frechet differ¬ 

entiable. In order to show this, we treat each of the four terms occurring in dsni) 
separately using the triangle inequality. In the following we sketch the calculation 
for the first two terms. In case of the first term, we have 


\m 






J Hp- 




ir2(s2) 


< 


47r 


I 2 —S — erjl 

\\V ^ 


H2(§2) 


We now use the definition of the ||•||^ 2 ^g 2 ^-norm given in (I48II on the last of these 
factors IIII in order to show that this term can be bounded by an ex¬ 
pression of the form 5(||(^||j:f2('s2) , ||? 7 ||j:f 2 (s 2 )) Il’ 7 llff 2 (s 2 )- Because the expression 
in (|48l) is quite long, we will also break it up in its individual summands and 
call them Ti,i for i = 1,...,3. We start with the first one sin 6d(pd6) for 

/ = ri^{(j),0) exp {—<1) — erj) (notice that we will drop the (tp, ^(-dependence from 
now on). In particular. 


Ti.i = 

< Ce 

< Ce 


(^j e sinOdipdO^ 


A 

ll-i-IU+dhll 


iL 


rf sin OdipdO ) 


ll-i-lloo+dhllo 


A / 

L4|| 2 _ ^g||</>|lo„+e||?j|| 


llir2(^S2) • 

This gives us a bound on the term Ti,i. In case of the second term, we have 
Ti 2 = j ^sin“^ 6 {dipf)^ + {def)^^ sinOdtfidO 

sin“^ 6 sm6d(pd6 


-L 


^1,2,1 


+ 


J (^2riderie ri‘^de{—(l){q) — €T]{q))e ^'n(Q)^ sinOdipdO 


^1,2,2 
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Using (a + 6)^ < 2 (a^ + b ^), we obtain that 

24 . 2,1 < J sm~^ 9 sin 6d(pd9 

<j sin“^ 0 ^|j277e“'^“''’’|| (d^rj)^ + — erj))^'^ dipdO 

^^IIJillL [ sm~^ 9 {d^r]f d(pde 

II lloo^A 

+ hill [ sm~^ 9 - en)f dipd9 

II " 007.4 

^Chlll hll?f2(s2) +<^11^1112(82) Ih + «’7lll2(s2) 

^<^hlll2(s2) +C'||??||^2(s2) Ih + e?7lll2(s2) • 

Similarly, we obtain for Ti. 2.2 that 

Ti,2,2 < C'||r?||^2(g2) Ih + '^^lll2(s2) • 

Using \^a + b < ^fd + \fb^ we obtain for ri .2 

+ 1.2 < c h||l 2 ( 82 ) (1 + U + ^v\U ■ 

Bounds for all other terms can be obtained similarly and will be omitted. □ 

C.2 Computing products of spherical harmonics 


Algorithm 2: Rules for the computation of products of spherical har¬ 
monics 


Require: a polynomial / of spherical harmonics of arbitrary degree. 
Ensure: an expansion in terms of spherical harmonics. 

1: If a is a constant, then define ^[a] := a and r[a] := a. 

2: Define Ni,^ := 

3: Define S [Yr] := Yr and T [Yr] := YT- 

4: If a is a constant, then define Sfay] := aS [y] and T[ay] ~ aT \y\. 
5: Define S [® -I- y] := S' \x\ -|- S [y]. 

6: Define S [aYT+p"] “ S [aS [YrY^^]]. 

7: If p > q and I — 2 > \m\, then define 


S [Uru;] := 


1 P Q ~l~ 1 ^l,m _ ^p,q o r-y^q 1 

l-m 2p + l A(p+i),, L P+i *-iJ 

_l_ 2Z — 1 P q Nl^rn _ ^p,q g 1 

^ I-m2p + l ^ 


Z — 1 + m Ni- 


I — m N(i_ 




-S [Yt^2Y,^] . 
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8: If I = m, then define 

S [YrVp^] := {-iyNi,^{2l - 1)!!S' [(1 - . 

9: If Z = m -I- 1, then define 

S [yryp"] := Ni,m{2l - l){-iy-\2l - S)!!^ [(1 - . 

10: If I — —m, then define 

S [yry;] := ^(2^ - 1)!!T [(1 - x^y/^Yy\ . 

11: If Z = —m -f 1, then define 

S [yry/] := - 3)!!r [(1 - x^)^xYy] . 

12: If I > \m\, then define 


S[a-xYr] ~ 


Ni, 


^(l + l),7n 

13: If I = \m\, then define 

S[a-xYr] := 

14: If / > m -I- 2, then define 

s [(1 - x^YYr] ■■= ^ 


I 771 “h 1 -^rn 


Nl,m o 


I 7TL -y^rn 


Ni, 


N, 




Ni, 


I -m + 1 
a —-T-:—y. 


2/ -b 1 


i-l-i 


21 + 1 

1 Nl^rn 


21 + 1 

15: If Z = m or Z = m -b 1, then define 

1 Nl^rn 

2Z -b 1 A^(( + l),(m,+ l) 


S [(1 - . 


S'[(1 - a;^)"yn :=- 


S [(1 - a:^)""^^^y(7i+^ . 


16: Define S [a(yr)"] := S [a(yr)''"^S [(yD^]] • 
17: If Z — 2 > Iml, then define 


S 


'(y” 


m\ 2 '] 21 — 11 — m+1 Ni . 


Ni, 


^ SlYr^YY^i] 


Z m 2Z -b 1 — 


2Z — 1 Z -b m f Ni^ 


—) s[y,!!iiy*!:!i] 

),m ) 


I —m 21 + 1 

^ 4 “ 1 ^l,m ^ 

-o [yi_2ri J . 


Z - m N^l- 2 ),m 

18: If Z = m, then define 

S [(yr)"] := iVi,^(-l)*(2Z - 1)!!S [(1 - xy^^Yr] . 

19: If Z = —m, then define 

S [(yr)"] - iV;,-(((2Z)!)"'(2Z - l)!!r [(1 - xy^^Yr] . 

20: If Z = m -b 1, then define 

S [{Yrfj := iVi,^(-l)*-i(2Z - 1)(2Z - 3)!!S [(1 - x^)^xYr] ■ 
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21; If I = —m + 1, then define 

5 [{YD^] 

■■= - 1)!)“^(2/ - 1)(2Z - 3)!!T [(1 - x^)‘~^xYr^ . 

22: If Z > |m|, then define 

S' [a • xYr] := 1+IHt j ^ 

23: If Z = |m|, then define 

o r ^l,m Z TTl “h 1 |. -tT-m i 

S [a • xYi ] := —^- T [aYi+j] . 

+ ^Z + 1 


24: If Z > —m + 2, then define 

(Z + m — 1)(Z + m)A'j 

,TitrT-Tr/-i 2\s —l/2-^m—11 

■ (2Z + 

25: If Z = —m + 1 or Z = —m, then define 


C.3 The bifurcation equation in terms of real spherical harmonics 


The bifurcation equation in terms of real spherical harmonics is given by 


2 

/real(M,A)(p) = ^ fm(u,\)Y2,m(p). 

m= — 2 
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, . 1 [ 1^2 1 (-l-hll207r)u_2«§ 

+ 448 - 12544(1 +32.) 


1S\U-2Uq 
49.(1 + 327r)2 


89 - 1792. + 250880.^) uhul 
1931776(1 + 32.)2 


1 98. + 4480. ) m_2«o (—1 + 1120.)-u_2'U-iMi 

25872(1 + 32.)2 ^ 6272(1 + 32.) 


36Au_2M-iWi V 2tv 1792. + 250880. ) m_iMi 
“ 49.(1 + 32.)2 965888(1 + 32.)2 

(23 — 9184. + 250880.^) M_2M_iWo'ni 
965888(1 + 32.)2 

5\/¥ (31 + 1120. + 50176.2) (_i + Ii20.)u2 

965888(1 + 32.)2 6272(1 + 32.) 

36Am2 2M2 (61 + 952. + 125440.2) u-2uiiU2 

^ 49.(1 + 32.)2 241472(1 + 32.)2 

y| (61 + 952. + 125440.2) u'i2UoU2 
120736(1 + 32.)2 ’ 

1 r -— (—1 + 1120.)rt— iUq 18Aw-iUo 

/_1 - M_1 + —V 12544(1 + 32.) ^ 49.(1 + 32.)2 

yi (577 + 5824.+ 1254400.2) 

5795328(1 + 32.)2 ^ V 

(—1 + 1120.)m?.iMi 36Au?.iMi 
6272(1 + 32.) 49.(1 + 32.)2 

y| (61 + 952. + 125440.2) uliuoui 
241472(1 + 32.)2 

yif (89 - 1792. + 250880.2) u-2ului 
965888(1 + 32.)2 

3\/¥ (111 + 672. + 250880.2) u-2U-iuf (_i + ll20.)M-2M-iri2 
965888(1 + 32.)2 “ 6272(1 + 32.) 

36Au-2rt-iW2 5y^ (31 + 1120. + 50176.2) u^_^U2 
49.(1 + 32.)2 + 965888(1 + 32.)2 


(221 + 12992. + 250880.2) m_2«_iMoW2 
965888(1 + 32.)2 

If (89 - 1792. + 250880.2) u'i2UiU2 
482944(1 + 32.)2 ’ 
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f , 1 2 (~1 + 11207r)MQ _ ISAttp 

/o - Mo + V ttmq ;^2544(1 + 32n) 497r(l + 32tt)^ 

(577 + 58247r + 12544007r2)Mo 1 _ 

5795 3 28(1 + 327r)2 

(—1 + 11207r)M_iMoMi 36 Am_iMoMi 
62 72(1 + 327r) 497r(l + 327r)2 

y| (577 + 58247r + 12544007r2) M-im^ui 
193 1 776(1 +327r)2 
yi (89 - 17927r + 25 0 8 807r^) jM? 

96 5 8 88(1 + 327r)2 

5^^ (31 + 11207r + 5017671^) U- 2 M 0 M? i _ 

965888(1 + 327r)2 

(—1 + 11207r)u_2MoM2 , 36 Am_2MoM 2 
62 72(1 + 3277) ^ 497r(l + 327r)2 

(31 + 112077 + 5017677^) m?.iMoM2 
965888(1 + 3277)2 


5^1 (31 + 112077 + 5017677^) M_2MgM2 
482944(1 + 3277)2 

(89 — 179277 + 25088077^) U- 2 U- 1 U 1 U 2 
965888(1 + 3277)2 

(-12365 - 205241677 - 1202073677^ + 104044953677®) m_2M-iMoMiM2 
199763277(1 + 3277)3(5 + 51277) 

(89 - 179277 + 25088077^) -u^aMi 
482944(1 + 3277)2 ’ 
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r \ 1 — (—1 -f 11207r)uorti 

h =Ami + — VSttmomi - , oo_^ 


ISAmqMi 


-VOTrtioWl-;-r--;-r^r 

48 12544(1 +327r) 497r(l + 327r)2 

577 -h 58 247r -b 12544007r2) ngwi ^ 11207r)M_iM? 

5795 3 28(1 + 327r)2 6272(1 -b 327r) 

36Au_iUi \f^ + 9527r -b 1254407r^) 

497r(l -b 327r)2 241472(1 -b 327r)2 


5^/^(31-bll207r-b 50 1 767r2)M_2iti i 


965888(1 -b 327r)2 224 V 2 

^ (89 - 17927r -b 2508 807r2) u-iulu2 (_i + ll207r)u-2MiM2 

965888(1 -b 327r)2 6272(1 -b 327r) 

36A-u_2Mi'U 2 -b 2508807r ) u_iMiU2 

497r(l -b 327r)2 965888(1 -b 327r)2 


^ (221 -b 129927r -b 2508807r^) m_2«oUiM 2 
965888(1 -b 327r)2 

'll (89 - 17927r -b 25088071^) u-2U-iul 
482944(1 -b 327r)2 ’ 


1 /l^ 2 V ^ (89 - 1792^ + 250880^2) 

/2 -Am 2 + —Wi 193 1 776(1 + 327r')2 


'448 V 2 ^ 193 1 776(1 + 327r)2 

g (89 - 179277-b 25 0 8 807r2)M-iu? i 

965888(1 -b 327r)2 Wa 


(—1 -b 11207r)7tQ7t2 18 AMoM 2 
12544(1 + 327r) 497r(l -b 327r)2 “ 

(—1 + 11207r)7t_lMlM2 36 Am_iMiM2 


V^UoU2 


I (l- 987r-b44 807r2)MgM2 

25872(1 -b 3277)2 


6272(1 -b 3277) 


4977(1 -b 3277)2 


23 — 918477 -b 25088077^) u-iUouiU2 
965888(1 -b 3277)2 


yH (61 -b 95277 -b 12544077^) U-2uiu2 (_1 + 112077)M_2Mi 
241472(1 -b 3277)2 6272(1 -b 3277) 

36A«_2«i (31 + 1120^ + 50176^2) 

^ 4977(1 -b 3277)2 + 965888(1 -b 3277)2 

yi (61 + 95277 + 12544077^) u-2Uoul 
120736(1 -b 3277)2 • 




60 


Michaela A. C. Vollmer 


C.4 Properties of the Bifurcation Equation 


Lemma 43. The bifurcation equation f{u,X) associated to the Euler-Lagrange 
equation for the Onsager free-energy functional is smooth. 


Proof. The bifurcation equation is defined by 


PTZ{u + v{u, A)) = f{u, A) 


where TZ is given by 


E{p,x) = /:{p,\) + n{p,\)- 


Since £ is a linear operator, the smoothness of E and TZ is the same. Moreover, 
the implicit function theorem tells us, that the smoothness of R and v is the same. 
Since the smoothness of the bifurcation equation is given by the minimal number 
of derivatives of TZ and v (using the chain rule) which are then equal, we conclude 
that the smoothness of / is the same as the smoothness of E (the Euler-Lagrange 
equation). The Euler-Lagrange operator is infinitely many times differentiable, see 
Proposition BU □ 


D Details of the dimension reduction based on invariant theory 

D.l The Wigner rotation matrix for degree I = 2 spherical harmonics 

The Wigner rotation matrix for complex spherical harmonics of degree / = 2 is 
given by 


/dll 

d2i 

dai 

d4i 

V dsi 


di2 di3 di4 di5 \ 
d22 d23 d24 d25 
d32 d33 d34 d35 
d42 d43 d44 d45 
ds2 ds3 ds4 dss / 


D{4>R, IpR, Or) 
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with 

dll 

di 2 = cos® sin , 

di 3 = cos^ (^) sin= (^) , 

di 4 = cos sin® , 

dis =e2'’^«-2'‘^«sin'‘ , 

d 2 i = _2e-‘^R-2»’^« cos® sin , 

d 22 = cos= (2 cos(6Ih) - 1), 

d 23 = cos cos{9r) sin , 

d 24 = e'’^«-''^«(2cos(0fl) -h l)sin= ('^) , 
d 25 = 2e®'’^«-'^« cos sin® j , 
d3i =y6e-2*’^«cos2(^)sin2(^), 
d 32 = —cos cos(6lfl) sin , 

d33 = (I COS^(6'i{) - , 

d34 = ■\/6e*’^” cos cos (dll) sin j , 
d35 =y6e2^’^«cos2(V)sin=(^), 

d4i = -2e^‘l>R-^i^R cos (^%-) sin® , 

d 42 = e'‘^«-"^«(2cos(0fl) -b l)sin= (%‘) , 
d 43 = —■\/6e*'^” cos cos(0_r) sin , 

d 44 = cos^ (%■) (2cos(6)fl) - 1), 

d 45 = cos® sin , 

dsi =e2'‘^«-2'’''«sin^(^), 
d52 = -2e2'‘^«-'’^« cos (^%-) sin® , 

d53=y6e2^^«cos2(^)sin2 (^) , 
d 54 = -2e^'‘^«+'’^« cos® sin , 

dss = cos'* . 

where 0^, ipR a-nd Or denote the Euler angles corresponding to a rotation. 


D.2 Representation of the group action for real spherical harmonics 

An application of the mapping given in (I31II yields an explicit form of the repre¬ 
sentation of the group action of 50(3) acting on R®, see (I34E 

/Mil Mi2 Mi3 Mi4 Mi5 \ 

M21 M22 M23 M24 M25 
M31 M32 M33 M34 M35 
A/41 A/42 A/43 A/44 A/45 
V M51 M52 M53 M54 M55 / 


M{4>r, tpR, Or) 
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with 


Mil 

Mi2 

Mi3 

Mi4 

Mis 

M21 

M22 

M23 

M24 

M25 

M31 

M32 

M33 

M34 

M35 

M41 

M42 

M43 

M44 

M45 

Msi 

M52 

M53 

M54 

M55 


cos{0n) cos(2i/)h) cos(2(^h) — i(cos(26ln) + 3) sin(2i/)fl) sin(2(^i{), 
sin(6»H)(cos(6»i{) sin(2V>fl) sin(0fl) - cos(2V>fl) cos(</)fi)), 
--y3sin^(6'n) sin(i/)H) cos(V>h), 

sin(6'H)(cos(6'H) sin(2i/'H) coa{(j)R) + cos(2i/)h) sin((?!>fl)), 

— cos(6Ih) cos(2i/)h) sin(2(^n) — i(cos(26ln) + 3) sin(2i^H) cos(20fl), 
sin(6li{)(cos(^/)H) 008 ( 2011 ) - 2 cos(6»i{) sin('0H) sin(0n) cos(0n)), 
cos(0n) cos(0h) cos(0h) — cos(20h) sin(0H) sin(0H), 

cos(0n) sin(0ii), 

— cos(6li{) cos(0n) sin(0n) — cos(26»n) sin(0i{) cos(0fl), 
sin(6'n)(- cos(0fl)) sin(20n) - i sin(20fl) sin(0H) cos(20i{), 
\/3sm^(6»n) sin(0n) cos(0h), 

\/3sin(0H) cos(6<h) sin(0H), 
i(3cos(26iH) + l), 

-ySsin)^!!) cos(0n) cos(0h), 
i\/3sin^(6»H) cos(20h), 

8111 ( 011 )(cos(efl) co8(0n) 8in(20ii) + 8in(0ii) cos(20ii)), 

CO8(20h) CO8(01i) 8in(0l{) + CO8(0 h) 8in(0ll) cos(0fl), 

—'y38in(0ll) CO8(0 h) CO8(0h) 

co8(20h) co8(0h) co8(0h) — cos(0ii) 8in(0ii) sin(0ii), 
i 8in(20i{) cos(0ii) cos(20ii) - 8in(0ii) sin(i/)n) 8in(20i{), 
|(cos(20ii) + 3) 008(201}) sin(20i}) + cos{9ii) sin(20i}) co8(20i}), 
sin(20i}) oos(20i}) sin(0i}) - sin(0i}) 8in(20i}) oos(0i}), 
iv^8in^(0i}) 008(201}), 

8in(0i}) sin(20i}) sin(0i}) - i 8in(20i}) co8(20i}) oo8(0i}), 
i(oos(20i}) + 3) 008(201}) oos(20i}) — cos { 6 r ) sin(20i}) sin(20i}). 


D.3 The explicit form of the Cartan representation 


The explioit form of the 5 x 5 matrix representing the Cartan representation is 
given by 


Me (01}, 01}, Or) 


/Cll C12 Cl3 Cl4 Cl5\ 
C2I C22 ^23 ^24 C25 
^31 C32 C33 Cii C35 
C4I C42 ^43 ^44 C45 
\ C51 C52 Cs 3 C34 C55 / 


with 


Cll = i (-4 008^(01}) 008^(01}) sin^(0i}) + (eos(20i}) + 2eos(20i}) - 1) eos(20i})) 
+2eos(0i})(cos(0i}) — 4eos(0i}) eos(0i}) sin(0i}) sin(0i})), 

C12 = 2eos(0i}) oos(0i}) sin(0i}) sin(0i}) eos(0i}) + cos(20r) eos^(0i}) sin^(0i}) 

+ (sin^(0i}) - 008^(01}) sin^(0i})) eos^(0i}), 

C13 = ^ ((cos(20r) + 3)eos(20i}) - 2sin^(0i})) sin(20i}) 

+ cos(0i}) cos(20i}) sin(20i}), 

Ci4 = 2oos(0i}) sin(0i})(sin(0i}) sin(0i}) - cos(0ii) cos(0i}) eos(0i})), 

C15 = -2eos(0i}) sin(0i})(eos(0i}) cos(0i}) sin(0i}) + cos(^r) sin(0i})), 

C21 = 008^(01}) sin^(0i}) sin^(0i}) + i cos(Or) sin(20i}) sin(20i}) 

(eos(20i}) + cos(20i}) (sin^(0ii) - cos'^(O r) oos( 20 i}))) , 
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C22 = - cos{(j)R) cos(V>fl) sin((/)H) sin(V'fl) cos^{O r) 

-l-|(cos(2(^fl) cos(2i/)fl) -I- 1) cos^{9r) 

-|-|(cos(26lfl) — 3) sin(2(?ifl) sin(2'^H) cos{9r) 

+ 5111 ^{Or) (cos^(</)h) cos( 2 V'fl) - sin^((^fl) sin^(i/>H)) , 

C '23 = -2cos{(f)R)cos'^{ipR)sm{<j>R)sm'^{0R) 

— cos( 0 h)(cos( 0 h) cos( 2 ?/)h) sin( 2 (?!>fl) -|- cos( 2 (/)h) sin( 2 i/)fl)), 

C24, = - 2 sin( 0 H) sin(i/>fl)(cos(V)fl)sin( 0 K) -|-cos( 6 'fl) cos( 0 fl) sin(i/>fl)), 

C25 = 2 8m{0R)sm{^R){cos{(j)R)cos{^R) - cos{9r) sm{(l>R) sm{ipR)), 

Csi = I ( - 4cos(6'fl) cos{2'4>r) sm{2(j)R) 

- {{cos{29r) + 3) cos(2(^h) - 6sin^(0H)) sin(2^/)fl)^, 

C 32 = I ^4cos(6lfl) cos(2V)fl) sin(2(?!)H) 

-I- ( 6 sin^( 6 lfl) -I- (cos(2^fl) -|-3) cos(2(^ij)) sin(2^/>fl)^, 

C33 = cos{0r) cos{2(j>R) cos{2iPr) — j{cos{20r) + 3) sin(20fl) sin(2^fl), 

C34 = sm{ 0 R){cos{ 2 ^R) sm{<j>R) + cos{ 6 r) cos{(j)R) sm{ 2 'ipR)), 

C33 = sin( 6 »fl)(cos( 6 lij) sin((?i_R) sin( 2 ^fl) - cos((^h) cos(2V>i?)), 

Cai = 5 sin( 0 H)(cos( 6 »i{)(cos( 2 (^fl) -|- 3) cos(V'fl) - 2 cos(<^fl) sm{(j>R)sm{iljR)), 
C42 = I sin(0H)(sin(20i{) sin(^/)fl) - cos{0R){cos{2(j>R) - 3) cos{'iPr)), 

C43 = sin{ 0 R){cos{ 0 R) cos{iPr) sm{ 2 (j>R) + cos( 20 h) sin(r^H)), 

C44 = cos{20r) cos{(f)R) cos(V’h) — cos{0r) sm{(l>R) sin(V>fl), 

C45 = cos{ 20 r) cos(V'fl) sin((/)i{) -|- cos{ 0 r) cos{(t>R) sm{ipR), 

Csi = sin(0H)(cos(i/'H) sin(20fl) -|- cos(0h)(cos(2(;!>h) -|- 3) sm{ipR)), 

C32 = I sin(0H)(cos(i/'H) sin( 2 (;/)_R) -|- cos( 0 h)(cos( 2 (;!)h) - 3) sin(V'H)), 

C33 = sin(6'fl)(cos(20H)cos(i/>fi) - 2 cos( 6 'fl) cos( 0 fl) sin( 0 H) sin(r/)i{)), 

C34 = — cos{0r) cos(V>i?) sin((^ij) — cos{20r) cos{cj)R) sm{iljR), 

C55 = cos{0r) cos{(f)R) cos{'iPr) — cos{20r) sm{(l>R) sm{ipR). 


D.4 Isomorphism between the representation of the group action of SO{3) 
acting on K® in (1341) and the Cartan representation 


Lemma 44 (Isomorphism between the group representations). An isomorphism 
between the two representations of the group action is given by 



0 

1 

0 

1 

\ 

0 

0 

1 

0 

-1 


1 

0 

0 

0 

0 


0 

0 

0 

1 

0 


Vo 

1 

0 

0 

0 

/ 


such that Me = ^. 

Proof. Omitted. 


(51) 


□ 


E Uniaxial solutions 


Lemma 45. The expression 


f / 

J 00=0 «/ ip 


277 

Vq=0 


K{0p,tpp,0q,ipq)p{0q)sin0q dOqdifq 


does not depend on pp. 
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Proof. Let R:, be a rotation around the z-axis. Using the rotational invariance of 
the interaction kernel k{p ■ q), we see that 


/ k{Rzp-q)exp{-(t){q))dq= k{p ■ Rlq) exp{-(t){q)) dq 
7s 2 y§2 

= k{p-q)exp{-(t){Rzq)) dq 
7§2 

= / ^)exp(-0(g)) dq. 

7§2 


Thus, the expression is invariant with respect to rotations around the z-axis and 
therefore independent of the variable ipp. □ 




